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ABSTRACT 


The  plane  wave  representation  is  generalized  in  terms  of  a  three-dimensional 
Fourier  transform,  to  yield  an  expression,  which  when  combined  with  the 
incident  field,  gives  the  total  electric  field  everywhere  inside  and  outside 
non- magnetic  scattering  bodies.  This  representation  requires  the  knowledge 
of  the  scattering  matrix  and  its  analytic  continuation  in  the  frequency  domain. 

The  use  of  short  pulse  data  to  determine  the  properties  of  uniformly 
coated  bodies  is  considered,  where  the  main  attention  is  applied  to  the  illuminated 
portion  of  smooth  convex  bodies  at  high  frequencies. 

Computational  results  are  obtained,  establishing  the  conditions  that  are 
required  to  determine  the  surface  of  a  conducting  body  from  knowledge  of  the 
total  near  field.  It  is  shown  for  the  cases  treated  (sphere  and  prolate  spheroid), 
that  two  separate  conditions  are  required.  One  of  these  conditions  yields  a 
single  surface  which  is  an  approximation  to  the  proper  surface,  whereas,  the 
other  condition  yields  a  set  of  surfaces,  among  which  is  the  exact  surface. 


iii 


TABLE  OF  CONTENTS 


ABSTRACT  iii 

I  INTRODUCTION.  1 

II  THE  PLANE  WAVE  EXPANSION  AND  ITS  APPLICATION  TO 

INVERSE  SCATTERING.  2 

ID  ASPECTS  ON  PULSE  SCATTERING  FOR  THE  DETERMINATION 

OF  UNIFORMLY  COATED  SCATTERING  BODIES.  17 

IV  ANALYSIS  OF  THE  COMPUTATIONAL  RESULTS  FOR  THE 
APPLICATION  OF  THE  BOUNDARY  CONDITIONS 

—  T  X  — T  =  0’  AND  |E.|  -  |Ej  =0.  55 

REFERENCES  88 

DISTRIBUTION  LIST 
DD  1473 


Iv 
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INTRODUCTION 

In  section  II,  the  plane  wave  representation  is  generalized  in  terms  of 
a  three-dimensional  Fourier  transform.  The  resulting  expression  when  combined 
with  the  incident  field  yields  the  total  field  everywhere  outside  and  inside 
non-magnetic  scattering  bodies.  This  representation  requires  the  knowledge 
of  the  scattering  matrix  and  its  analytic  continuation  in  the  frequency  domain. 

The  main  emphasis  is  placed  upon  homogeneous  bodies,  a  special  use  of  which 
is  the  conducting  body,  and  various  relations  are  obtained.  The  extension  to 
non-homogeneous  bodies  is  briefly  considered.  The  results  presented  here  are 
not  complete,  since  a  significant  amount  of  analysis  has  yet  to  be  performed. 

The  possibility  of  correlating  monostatic  short  pulse  data  and  bistatic  data 
for  a  fixed  frequency  and  direction  of  incidence  may  be  achievable. 

The  use  of  short  pulse  data  to  determine  the  properties  of  uniformly  coated 
conducting  bodies  is  considered  in  section  III.  The  main  attention  is  focused 
upon  the  illuminated  portion  of  smooth  convex  bodies  in  the  high  frequency  case. 

Finally  in  section  IV,  computational  results  are  presented.  It  is  shown 
for  the  two  cases  of  conducting  surfaces,  the  sphere  and  prolate  spheroid, 
that  the  two  conditions  |e  S|  -  |  E  1 1  =0,  and  E  x  E*  -  0  are  required  to 
establish  the  surface  of  the  body  from  knowledge  of  the  total  near  field.  The 
first  condition  yields  a  single  surface  which  is  an  approximation  to  the  proper 
surface.  The  second  condition  yields  a  set  of  surfaces,  among  which  is  the 
proper  surface. 
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II 


THE  PLANE  WAVE  EXPANSION  AND  ITS 
APPLICATION  TO  INVERSE  SCATTERING 


It  has  been  shown  previously,  that  if  the  scattered  far  field  (phase, 

amplitude  and  polarization)  is  known  over  all  directions  of  observation,  for 

a  fixed  direction  of  incidence,  then  a  particular  class  of  plane  wave  representations 

yielded  expressions  for  the  near  field  which  held  down  to  the  surface  of  smooth 

convex  bodies.  The  plane  wave  representations  involved,  contained  the  far 

scattered  field  as  an  analytic  function  of  the  angles  of  observation.  The  use  of 

the  plane  wave  representation  will  be  generalized  here,  to  yield  expressions 

that  hold  everywhere  in  space  including  the  interior  of  the  scattering  body. 

To  begin,  the  analysis  will  be  restricted  to  non-magnetic  bodies  (although 

it  could  be  easily  generalized  to  include  such  cases)  and  the  geometry  of  the 

scattering  body  will  be  limited  by  placing  certain  restrictive  analytical  properties 

on  the  surface  S  which  encloses  the  volume  V  of  the  scattering  body.  The 

3 

classes  of  surfaces  chosen,  will  belong  to  class  C  defined  by  Barrar  and 
Dolph  as  follows: 

3 

A  surface  S  is  said  to  belong  to  class  C  if  there  exists  a  finite  number 

m  of  images,  x'  =  x^u,  v),  yV  =  y^u,  v),  z  -  zy(u,  v),  v  =  1.2.3. . .  m  of  the 
2  2 

disk  u  +  v  <  1  that  cover  the  surface  S  and  such  that  the  third  derivatives  of 
x' ,  yV  and  z  with  respect  to  u  and  v  exist  and  are  continuous. 

Harmonic  time  dependence  exp(-iut)  will  be  taken  in  which  case  Maxwell's 
equations  become 


e  V  a  H 
o  —  — 


-i(E/k)  ' 


k  outside  the  body 


where  k  ,  the  square  of  the  propagation  constant  in  the  body,  is  associated  with 


2 


2  2 

the  free  space  propagation  constant  k,  i.e.  =k^(k). 

The  incident  electric  intensity  will  be  expressed  in  the  form 


i  ,  v-3/2  ik  •  x 

E  =  (27t)  '  a  e - 


(2.1) 


where  k  is  the  direction  of  the  incident  wave,  and  a  ,  the  unit  vector  denoting 

A 

polarization,  will  depend  upon  both  k  and  index  j  (i.e.  a  =  a  (k,  j)  where  the 

* 

index  j  takes  on  the  value  1  if  the  polarization  is  in  the  0  direction  of  a  fixed 

A 

coordinate  system,  and  takes  on  the  value  2  if  the  polarization  is  in  the  ^ 
direction. 


Homogeneous  Body 

The  homogeneous  body  will  be  treated  first.  In  this  case,  it  follows 
from  Banar  and  Dolph,  that  the  total  electric  intensity  E  induced  by  the  plane 
wave  Eq.(2. 1),  incident  upon  the  body,  satisfies  the  following  integral  equation 
(which  has  a  unique  solution) 


E(x) 


0(E  •  n  )  ds 


(2.2) 


where 


ikR 

0  «  ,  RMx-x') 


and  E  is  the  value  of  E  obtained  by  approaching  the  surface  from  the  interior 
of  the  body.  If  E+  is  the  value  obtained  by  approaching  the  surface  from  the  exterior 
of  the  body,  it  follows  that 


2  2  + 
kx  E  •  n  =  k  E 


(2.3) 
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In  addition,  the  magnetic  field  can  be  expressed  in  the  form 


H(x)  =  H1(x)  + 


(krk2) 


—  —  47riWM 


V  A  /  E  0  dv  . 


(2.4) 


Before  deriving  a  plane  wave  representation  for  E(x),  the  following 
vector  will  be  introduced 


T(k',k)  = 


CM 


(2tt) 


3/2 


—  i  k'  *  x'  ,  ...  .  ik'  /  -ik'-x'  -  ,  . 

e  —  ~  E (x  ,  k)  dv - /  e  —  —  E  •  n  ds 

^  V  "S  (2.5) 


From  Eq.(2.2),  it  is  seen  that  when  |x|  — ►  oo  ,  in  the  direction  given  by 
the  unit  vector  k' ,  then  the  total  field  becomes 


iklxl 


E  (x)  ^ 

I  x|  -»00 


-i(-)+/?liT  x-'- 


(2.6) 


indicating  that  the  vector  T  (k\  k)  is  related  to  the  scattering  matrix.  When 
|k'|  =  | k|  ,  T(k',k)  is  a  measurable  function  being  proportioned  to  the  far  scattered 

field  in  direction  k'  ,  produced  by  a  plane  wave  of  frequency  ck,  incident  upon 

A 

the  body  with  direction  of  incidence  given  by  k  .  The  following  theorem  may 
now  be  proven. 

Theorem: 


If 


t  (k',k)  =  — - 


(27 r) 


3/2 


-ik'*  x  „  ,  , .  , 

e  -  -  E (x,  k)  dx  , 


I(k',k) 


then  ^(k1  •  k)  =  6  (k  -  k1)  a  -  lim  — - - „ - 

e-*  o  k  -  k’  +  ie 


(2.7) 


(2.8) 
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Proof: 


It  follows  from  Eq.(2.2)  that 

(krk2) 


6  (k',  k)  =  6  (k-k')  a  +■ 


477 


(2  it) 


3/2 


E  (x',k)  I  e  1~  —  0  dx  dv 


\  1  6 


-ik'-  x 


(27 t) 


3/2  - 


V  •  /  0  (E*- n)  ds  >  . 


On  setting 


\jj  -  J  0(E  •  n)  ds 
S 


one  can  show  that 


-ik'  -  x 

e  —  —  V  ^  dx  =  lim 


-ik'-  x  ,  . ,  ,  /  -ik'-  x  . 

e  —  —  ^/ds+ik'  e  ~  —  ip  dx 


R  — MX)  w 

oo  v  S 


00 


when  S  is  the  surface  of  a  sphere  of  radius  R  .  Letting  k  have  a  small 
oo  oo 

imaginary  part,  it  is  seen  that  the  surface  integral  will  vanish  when  R^  oo  . 
The  resulting  integrals  can  be  reduced  as  follows 


ei]D  •  (x  •  x')  -  ik'  -  x 


2 

P 


-  k  -  ie 


-ip  •  x' 
e 


2  i2  • 

p  -  k  -  ie 


6  (j>  -  k') 
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2  2 
k  -k'  +ie 


e 


Combining  the  above  expressions,  one  obtains  the  result 


t_  (k'-  k)  =  6  (k-k')  a  -  lim 

6  — ►  o  2  2 


k2-k'2  +  ie 


An  immediate  consequence  of  the  above  theorem,  is  the  following  corollary 
which  expresses  the  total  field  everywhere  in  space  in  terms  of  the  incident 
field  and  the  quantity  T  (k',  k)  , 


(2.9) 


If  the  integration  space  ja  is  expressed  in  spherical  polar  coordinates  (p,  0^,  0^) 
where  the  range  of  the  variables  are  -co<p<oo,  0  <  6^<  7t/2  ,  and 
0  <  0p£  27r  ,  then  in  the  above  integral,  the  contour  of  the  variable  p  bends 
above  the  pole  at  p  =  -k,  and  below  the  pole  at  p  =  k  . 

It  can  be  shown  in  a  similar  manner  that  the  magnetic  field  can  be 
expressed  in  a  similar  form 


(2.10) 


Integral  Equation  for  T  (k',  k) 

Applying  the  condition  that  V  •  E  =  0  in  the  volume,  one  can  express 
Eq.  (2.5)  in  the  form 


6 


dx 


(2.11) 


where  the  surface  integral  is  replaced  by  a  volume  integral.  On  substitution  of 
Expression  (2. 9)  for  E,  Eq.  (2.11)  reduces  to  the  integral  equation 


Defining  V(k',k)  to  be  the  Fourier  transform  of  the  body,  as  follows 


V(k\k) 


(2  7T  )3 


i(k'-k)  •  x  . 
e  —  —  —  dx 


(2.12) 


one  obtains  the  following  fundamental  integral  equation  for  T  , 

\ 


(M 


a  V (k,  k')  + 


T(£,  k) 

p  -  k 


-  T(k\k)  - 


k'  (k'  •  T) 

2  2 
k^  -  k 


(2.13) 


For  the  case  of  almost  transparent  bodies,  i.e.  where 
2  , 2' 


(k2  -  k2)  V  (k,  k') 


is  small,  an  approximate  solution  to  the  integral  equation  is  obtained  through 
a  process  of  integration,  yielding  the  Born  approximation 


7 


T  (k',k)~A(k',k)  * 


aV(k,k')  +  /dp 


whereA(k',k)  is  the  operator 


A  <£’  k)  V  ( k,  p)  a  +  .  .  . 

p  -  k 


A 


1  - 


k'(k'-  ) 


For  the  case  of  a  perfect  conductor,  i.e.  .  Imk^->oo ,  the  integral  equation  (2. 13), 
cannot  be  used  by  itself.  An  additional  equation  is  required.  The  details  on 
the  particular  aspect  are  given  later  on. 

Relationships  Among  T  (k1,  k) 

In  order  to  attempt  to  reduce  the  number  of  measurements  needed  to 
determine  the  scattering  body,  or  in  order  to  correlate  different  sets  of 
measurements  (such  as  for  different  directions  of  incidence),  various  relationships 
involving  T  (k  \  k)  for  different  k'  and  k  are  sought. 

Denoting  the  complex  conjugate  of  a  quantity  be  a  super  bar,  it  can  be 
shown  that 


V(k,  k')  =  V  (-k,  -  k') 


Hence  it  follows  that  the  complex  conjugate  of  Eq.  (2. 13)  yield  for  k^  real 


CM  * 


V(-k,-k')  +  /  d£  1<£.W 


T(k'  k)-JK!S12> 

k»  -  k 


P  -  k 

'  (2.14) 

In  the  above  equation,  replace  k  by  -k  ,  and  k'  by  -k'  .  In  doing  so,  the 
polarization  vector  a  must  be  changed  appropriately,  i.e. 


a  (k,  j)  =  a  (-k,  j) 
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where 


=  1  for  j  =  1 
=  -1  for  j  =  2 


This  results  in  the  following 

(k^-kV  e  a  (-k,  j)  V  (k,  k')  +  /  d£  T(-£,-k)  - 

P  -  k 

k'  -  k 


(2.15) 


It  follows  that  on  subtracting  Eq.  (2. 15)  from  Eq.  (2. 14)  for  j  =  1,  and  adding 
them  for  j  =  2,  that  the  vector 

T  (k\  k,  j)  -  c  T  ( -  k ' ,  - k,  j) 

is  a  solution  of  the  homogeneous  equation  corresponding  to  Eq.  (2. 13). 

Since  the  solution  is  unique,  the  following  must  hold 

-  2 

T(k\k,  j)  -  e  T(-k',-k,  j)  ;  kj  real  (2.16) 

2 

This  result  holds  only  for  k  ,  real  and  finite.  However  it  can  be  shown  that 

it  will  also  apply  to  the  case  of  a  perfect  conductor,  where  Im  k  ->  oo  . 

2  ^ 

For  general  real  or  complex  values  of  k^  ,  an  additional  set  of  relations 
may  be  obtained.  To  proceed,  we  first  note  that 


V(£,k')  V(k\k) 


V(jD,k)  , 


9 


hence  upon  operating  on  Eq.  (2.13)  with  respect  to 


dk'  V  (k\q) 


it  is  seen  that  the  left-hand  side  is  invariant.  Defining 
T(k\k)  = 

one  immediately  obtains  the  result 

/  d£V(£,k')  TCp.k)  -  T(k\k)  .  (2.18) 


k(k*-  ) 

2  2 
k'  -  k 


T  (k',k) 


(2.17) 


This  equation  will  be  shown  to  the  additional  equation  required  for  the  case  of 
a  perfect  conductor.  For  present  purposes,  this  equation,  combined  with 
Eq.  (2. 14)  which  can  be  expressed  in  the  form 

r  \ 


(k2  -  k2) 


a  V  (k,  k’)  + 


V(f,~')  I  (P.  » 

p  -  kz 


•  =  T(k',k) 


(2.19) 


can  be  used  to  eliminate  the  volume  factor  V(k',  k)  .  In  Eq.  (2. 19)  replace 
k'  by  -g  ,  and  operate  on  the  equation  with  an  arbitrary  unit  vector  b  ,  yielding 


b-a  V  (g,  -k) 


v  (a.  -JB> 

2  ,2 

p  -  k 


b  •  T(jg,k) 


b  ■  T(-g,k) 


Operate  on  thiB  equation  with  the  integral  operator 

r  ~ 

/  dq  T(g,k') 


10 


and  employ  relation  (2.18)  to  yield 


(2.20) 


Since  b  is  an  arbitrary  unit  vector,  this  equation  constitutes  three  independent 
equations.  For  the  case  of  a  perfect  conductor,  the  right-hand  side  of  Eq.  (2.20) 
vanishes.  It  can  be  shown  that  for  the  perfectly  conducting  case,  Eq.  (2.20)  will 
yield  the  law  of  reciprocity  when  |k'|  =  jk|  . 

Extension  of  the  Integral  Equation  to  the  Perfectly  Conducting  Case. 

2 

Before  considering  the  case  where  Im  k^-*ao  ,  an  analysis  upon  the 
decomposition  of  the  field  into  the  interior  and  exterior  components  will  be 
treated.  The  vector  T  will  be  decomposed  into  parts  T+  and  T  as  follows: 


E*  + 


d  2  = 


'  E  (x,  k) 
0 


for  x  outside  V 
for  x  inside  V  , 


T  dja  = 


E  (x.k) 


0 

\ 


for  x  inside  V 
for  x  outside  V 


It  then  follows  from  Eq.  (2.21)  by  integrating  over  V  ,  that 


(2.21) 


(2.22) 


a  V(k,k')  + 


V  (jd,  k') 


T+  (jg,  k)  dja  =  0 


(2.23) 
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and  from  Eq.  (2.22) 


T“  (k\k) 

— 2 - T 

k'  -  k 


(2ti) 


3/2 


E  (x,  k)  e  —  dx 


=  aV(k,k')+  /  d_p  V2^^r  =  ./ d-E  ~4M=- I  k) 

p  -  k 


"2  -  ,T 

p  -  k 


Applying  these  results  to  Eq.  (2. 13),  one  obtains 


(2.24) 


k'2  _  k^  \  /v/ 

1  (k'.is)  =(— - 2~J  I  (k'.k) 


(2.25) 


+ 


from  which  the  following  relationship  between  T  and  T  may  be  derived 


T"  (k'.k)  = 


V-k2 

k2  -  k'2. 


+ 


<  T  (k'.k)  - 


(k2-k2) 


+ 


2  a  .2  .2^  -  -  - 


k"  (k-2-  kO 


k'  (k'  •  T  )>  .  (2.26) 


Hence  on  eliminating  T  from  Eq.  (2.24),  one  obtains  the  additional  equation 


for  T 


T+  (k’,k)- 


(A  -  V 


k1  (k1  •  T+ ) 

k’2  -  k2 


djD  V(jD.k') 


G-2) 


X*{£,k)-  y  2  2  £(£-l+» 

k2  <p2-  k2) 


(2.27) 
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The  case  of  perfect  conductivity  may  be  considered  by  letting  Im  -*oo  . 

It  is  seen  from  Eq.  (2.25)  that 

T"  (k',k)  - 

i.e.  T  ~  T+  + 

Eqs.  (2. 23)  and  (2.27)  reduce  to  the  form 

(fi,k)  =  0  (2.28) 

T(k',k)  =  J  d£  V(j3,  k’)  T(j3,k)  .  (2.29) 

These  two  equations  constitute  the  two  integral  equations  that  are  required  for  the 
case  of  the  perfect  conducting  body.  Eq.  (2.28)  insures  that  the  total  fields 
vanish  inside  the  body,  whereas  it  can  be  shown  that  Eq.  (2.29)  corresponds 
to  Maue’s  integral  equation.  The  problem  with  working  with  Maue's  integral 
equation  by  itself,  is  that  it  does  not  always  yield  a  unique  solution,  since  there 
exists  a  set  of  resonant  frequencies  such  that  non- vanishing  interior  fields 
satisfy  the  condition  n  x  E  vanish  on  the  surface  of  the  body.  Thus  the  additional 
Eq.  (2.28),  removes  this  uniqueness  problem  by  requiring  that  the  field  vanish 
inside  the  body. 

Alternative  Representation  Exterior  to  the  Body 

For  a  fixed  direction  of  incidence,  the  quantity  T(k',  k)  is  required 
for  all  values  of  k'  ,  in  order  to  obtain  expressions  for  the  field  everywhere  in 


a  V  (k,  k')  + 


2  - 
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space.  However  it  will  be  shown  that  restricting  the  requirements  on  the  knowledge 

of  T(k',k)  to  the  values  of  k'  which  lie  on  the  sphere  |  k'|  =  |k|  ,  the  total 

field  can  be  obtained  everywhere  outside  the  minimum  convex  surface  enclosing 

the  body.  To  show  this,  let  the  plane  z  =  ZQ  be  the  tangent  plane  of  the  body, 

such  that  the  body  lies  in  the  half- space  z  <  zq  .  For  points  in  the  half- space 

z  >  z  ,  the  following  representation  of  the  near  scattered  field 
o 


ES(x) 


(2/)3/2 


T(£,k)  d  2 


(2.30) 


can  be  reduced  as  follows: 


Set  p  =  k  sin  a  cos  8 
x 


p  =  k  sin  a  sin  8 

y 


pz  --  kq 


when  the  domains  of  integration  are 

0  <  8  <  2i r  ,  0<Q'<7r/2-ioo,  and  -  oo  <  q  <  oo 

Expression  (2. 30)  becomes 


E 


(2n) 


3/2 


g(&,/3,  q)  cosQ'Sino'dQ'd]3d( 
(2.31) 


From  expression  (2.5)  ,  it  can  be  seen  that  for  z  >  zq  ,  the  contour  in  the 
following  integral 
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e  q  T  (p,  k)  dq 

2  2 
q  -  cos  a  -  ie 


may  be  deformed,  to  yield  the  following 
ik  cos  c?  z 

7r  i  e  T  (k  sin  a  cos  ]3,  k  sin  a  sin  /3,  k  cos  a  ;  k) 

Hence  expression  (2.31)  may  be  placed  in  the  form 
7r / 2— i  oo  /  27r 


E 


ik 
2  7T 


e*~  Eq (ch,  j8)  sin  a  d  a  d  P 


(2. 


where 


k'  =  k  (sin  a  cos  jS,  sin  a  sin  cos  a) 


and  E  is  the  amplitude  and  phase  of  the  far  scattered  field  in  direction 
given  by  (0,  0),  i.e. 

ikR 


Expression  (2.32)  was  developed  by  alternative  techniques  in  earlier  reports 
(Weston,  Bowman,  1966). 

The  Inhomogeneous  Body 

The  results  may  be  extended  to  include  inhomogeneous  non-conducting 
2  2  2 

bodies,  i.e.  where  k^  =  u  eiu  (a  real  function)  and  k^  varies  continuously 
in  the  medium.  In  this  case,  the  appropriate  integral  equation  for  the  total 
electric  field  is  (Barrar  and  Dolph) 


32) 
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which  has  been  shown  to  possess  a  unique  solution,  although  its  existance  has 
not  yet  been  proven. 


For  the  homogeneous  body,  it  can  be  shown  that 


In  this  case,  the  appropriate  integral  equation  becomes 


~  2  #  2  /  V  (n  k'l 

T  (k\  k)  =  k  aV  (k,k’)  +  k  /  T(jD,k)dj3 

P  -  k 


where 


V*  (k,  k')  =  — 

(  27T  ) 


dv 


with  k^  =  k^  N  (x)  . 
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ASPECTS  ON  PULSE  SCATTERING  FOR  THE  DETERMINATION 
OF  UNIFORMLY  COATED  SCATTERING  BODIES. 


3. 1  Introduction. 

A  convenient  method  for  the  identification  of  uniformly  coated  scattering 
bodies  may  employ  the  application  of  pulse  scattering  since  more  information 
may  be  obtained  by  such  a  transient  analysis  than  from  steady  state  results. 

The  distortion  of  a  pulse  returned  from  a  scattering  obstacle,  in  fact,  does 
yield  some  additional  information  pertaining  to  the  properties  of  the  scat¬ 
tered  e.g.  the  behavior  of  the  leading  wavefront  of  a  scattered  pulse  usually 
indicates  something  about  the  composition  of  the  body,  whereas  the  behavior  of 
the  trailing  return  of  a  scattered  pulse  is  related  to  the  shape  of  the  body  and 
its  radii  of  curvature. 

In  this  preliminary  treatment  on  the  identification  of  coated  scatterers, 
the  merits  and  the  restrictions  of  such  a  transient  approach  will  be  investigated, 
using  the  solution  of  pulse  scattering  by  planar  layered  structures  where  both 
the  mono-  and  bi- static  solutions  of  the  scattered  field  will  be  considered. 

Here  the  wavefront  technique  is  employed  to  solve  the  problem  of 
scattering  of  a  finite  rectangular  pulse  from  a  lossy  dielectric  slab  mounted 
on  a  perfectly  conducting  planar  surface.  The  general  case  of  oblique  incidence 
is  treated,  where  in  particular  the  behavior  at  normal  and  at  critical  (Brewster's 
angle)  incidence  are  demonstrated. 
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From  the  presented  theoretical  results  practical  methods  will  be  deduced 
which  may  be  applied  to  the  problem  of  identification  of  the  properties  of 
uniformly  coated  scatterers. 

3.2  Theory 

3.2.1  The  Laplace  Transform  of  Maxwell's  Equations  and  the  Boundary 
Conditions. 

The  solution  of  a  non-steady  state  problems  in  electromagnetic  theory  is 
facilitated  by  a  transform  type  approach.  Here  the  Laplace  transform  will 
be  used  since  it  is  defined  for  a  wider  class  of  functions  than  the  Fourier 
transform,  in  particular  to  the  extension  to  more  complicated  source  dependence. 

The  Laplace  transform  of  a  function  f(t)  is  defined  by  the  relation 


f00 

f(t)  e  St  dt 

•  0 


F(s) 


(3.2.1) 


The  inverse  transform  from  the  s  domain  to  the  t  domain  is  accomplished 
by  means  of  the  relation 


•B  +i  oo 

F(s)  eS*  ds 


JB  -ioo 


(3.2.2) 


where  B,  a  positive  constant,  is  chosen  such  that  F(s)  is  analytic  for 
Re  I  a  I  >  B.  In  particular  the  relation  holds: 


n-1 

s°  F  ( s)  -  ^ 

fco 


sn"i-1  F%) 


(3.2.3) 


The  vector  wave  equation,  satisfied  by  the  transformed  field  vector  _  (r,  s) 
is  given  by  2 

V  x  V  x  t  (r,  s)  +  ~  (1  +  — )  ^(r,  s)  (1  +  — )  E  (r,  0) 

—  —  2  es—  2  es  —  — 

c  c 


2  at 


c 


E  (r,  0) 


(3.2.4) 


where  ne  -  ijc' 
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The  equation  of  the  transformed  magnetic  field  vector  (r,  s)  is 
obtained  by  replacing  by  £  and  E  by  H.  Since  the  quantity  of  interest 
here  is  the  scattered  field,  it  is  always  possible  to  define  the  time  origin 
such  that  the  initial  values  of  the  scattered  field  are  zero  everywhere  in 
space,  and  thus  the  transform  of  the  scattered  field  vector  satisfies  the 
homogeneous  vector  wave  equation: 


V  x 


(3.2.5) 


The  boundary  conditions  satisfied  by  the  em  field  vectors  at  the  interface  of 
two  media  (the  normal  to  the  interface  n,  a  unit  vector,  is  directed  from 
medium  (1)  to  (2))  are  given  by  the  familiar  relations: 

a)  Both  media  of  finite  conductivity: 

n  x  j^E2^£»  t)  -  E j(r,  t)  j  =  0  or  n  (r,  s)  -  ^(r,  s)j=  0 

n  x  [h2^-’  =  0  or  n  x  s)  -  2^(r,  s)  =  0  .  (3.2.6) 

b)  Medium  (1)  being  a  perfect  conductor: 

n  x  (r,  t)  :  0  or  n  x  8)  =  0 

n  x  Hg (r,  t)  =  K (r,  t)  or  n  x  ^  (r»  8)  =  -^(r,  s)  .  (3.2.7a) 


In  cases  where  the  approximated  Leontovich  boundary  conditions 


-n  x  (  n  x  E  (r,  t)  =  Z  n  x  H  (r,  t)  (3.2.7b) 

with  Z,  the  surface  impedance,  may  be  implied,  the  transient  problem  must 
be  handled  with  care,  since  the  error  introduced  in  the  solution  of  the  diffraction 
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U  8  v 

of  a  harmonic  wave  (k  =  —  =  i—  )  by  a  curved  obstacle  of  curvature 

p  and  refractive  index  n,„  is  of  the  order  of  (n,n^)  .  Thus  only  the 

12  12  c 

behavior  of  the  scattered  fields  in  the  vicinity  of  the  wavefronts  which  is 
related  to  the  asymptotic  behavior  of  the  transformed  field  quantities  for 
large  values  of  the  transform  variable  s,  can  be  properly  derived.  Con¬ 
sequently  the  application  of  the  Leontovich  type  boundary  condtions  will  lead 
to  erroneous  results  for  the  field  vectors  when  the  time  delayed  from  the 
arrival  of  the  wavefront  is  large. 

3.2.2  Expansion  of  the  Transformed  Field  Quantities  as  a  Series  of 
Time -delayed  Terms  Determined  by  the  Optical  Wavefronts. 

In  Steady  state  electromagnetic  theory  it  has  been  shown  that  the 
field  scattered  from  an  obstacle  can  be  obtained  by  summing  the  contribution 
due  to  the  various  optical  rays  associated  with  the  scatter  geometry  of  the 
obstacle.  In  transient  analysis  such  families  of  rays  may  in  addition  be 
characterized  by  time  delay  factors  due  to  the  geometry  of  the  obstacle. 

To  verify  this  aspect  two  different  types  of  problems  will  briefly  be  con¬ 
sidered. 

3.2. 2.1  Pulse  Scattering  from  a  Finite  Obstacle 

Consider  the  configuration  of  Fig.3-1  where  a  plane  wave  is  incident 
on  a  finite  obstacle  of  closed  surface. 


FIG.3-1:  GENERALIZED  RAY  SYSTEM. 
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Q 

The  scattered  ray  E  Is  the  conventional  geometrical  optics  result,  where  the 


position  of  the  specular  point  T  is  determined  by  the  requirement  that  the 
optical  distance  between  reference  plane-scatterer-observation  point  is  a 


minimum.  In  addition  two  families 


the 


creeping  wave  contributions,  must  be  considered,  where  the  mth  term 


associated  with  the  ray  which  undergoes  m  complete 


in  the  sum 


circulations  around  the  obstacle  and  then  travels  to  the  observation  point. 

The  path  on  the  obstacle  taken  by  each  of  these  rays  is  determined  by  the 

condition  that  the  optical  distance  be  a  minimum. 

Since  wavefronts  of  non -monochromatic  waves  satisfy  geometrical 

optics,  the  above  described  diffraction  process  in  conjunction  leads  to  the 

hypothesis  that  the  solution  of  the  pulse  diffraction  problem  will  be  facilitated 

if  the  terms  corresponding  to  the  various  monochromatic  optical  rays  are 

Identified  and  treated  separately.  The  result  is  expressed  in  the  form  of 

an  infinite  sum  of  residue  and  branch  cut  contributions,  where  the  individual 

oo 


contributions  must  be  summed  in  a  manner  such  that  the  rays 
are  absent  for  t  <  T  . 


m=n 


n 


3 . 2 .  2 .  2  Pulse  Scattering  from  Layered  Structures. 

Another  class  of  problem  is  demonstrated  in  Fig.  3-2,  the  problem  of 
pulse  scattering  from  a  nonperfectly  conducting  slab.  The  incident  pulse  will 


TT 


i 


Ref.  Plane  of 
incident  wave 


E 


s 


3 


FIG.  3-2:  PULSE  SCATTERING  FROM  A  DIELECTRIC  SLAB. 
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be  decomposed  into  an  Infinite  sum  of  partial  pulse  returns,  where  each 
of  the  partial  returns  are  delayed  by  a  time  factor  T  ,  depending  upon  the 
material  properties  and  the  depth  of  the  layer.  The  amplitudes  of  individual 
partial  returns  are  obtained  by  selecting  the  proper  coefficient  of  the  expanded 
CW  reflection  coefficient. 

3. 2.  2.  3  Geometrical  Optics  Interpretation. 

From  the  two  presented  examples  it  consequently  seems  that  the  solution 
of  pulse  diffraction  problems  can  be  handled  most  efficiently  by  treating  each 
term  in  a  wavefront  expansion  separately,  where  the  amplitudes  of  individual 
wavefronts  are  given  by  the  corresponding  term  of  the  expanded  CW  scattering 
coefficient. 

The  solution  of  the  vector  wave  equation  (2.5)  for  the  transient  case 
is  obtained  in  the  same  manner  as  the  solution  of  the  vector  wave  equation 
for  monochromatic  fields.  In  general,  the  transform  of  the  field  vector 
^8(r,  s)  can  be  expressed  as  the  result  of  the  same  vector  operations  on  two 
scalar  functions  which  are  solutions  of  a  scalar  wave  equation  Consequently, 
in  discussing  the  decomposition  of  the  field  into  its  optical  components,  it 
will  suffice  to  investigate  the  behavior  of  a  scalar  function  F(r,  s)  which 
satisfies  the  scalar  wave  equation: 

2  s2  a 

V  F(r,  s)  -  -r-  (1  +  —  )  F  (r,  s)  =  0  (3.2.8) 

—  Z  es  — 

c 


where  F(r,  s)  will  be  rewritten  into  the  form  of  a  sum  of  terms  which  can 
be  identified  with  the  partially  delayed  optical  wavefronts: 

oo  oo 


F(r,  .)-  r  F  (r,  a)  =  ZI!  8>  e 

1 — rf  m  6 — rf  m 


-Sd  +—  )1/2T°  (r) 
es  m  — 


m=0 


m=0 


(3.2.9) 


Thus  the  following  relation  is  obtained: 
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Operating  on  F  (r,  s)  and  equating  the  result  to  zero,  it  is  found  that  the 


m 


functions  U  (r,  s)  and  T  (r)  must  satisfy  the  relation: 
m  —  m  — 


2 .  or . 
s  (1  +  — ) 
es 


VT°  (r) 
m  — 


2  J_ 

2 

c  _l 


U  (r,  s)  -  s  (1  +  —  )' 
m  —  es 


:i  +  — >1/2 


2VU  (r.  s) 
m  — 


VT°  (r)  +  U  (r,  s)  V2  T°  (r) 
m  —  m  —  m  — 


+  V  U  (r,  s)  =  0  (3.2.11) 

m  - 


where  the  spatial  time  factor  T  (r)  must  satisfy  the  Eikonal  equation 


V  T  (r) 
m  ~ 


-  Uc 


(3.2.12) 


If  s  is  replaced  by  s  =  -iu  and  or  =  0,  the  remainder  of  Eq.(3 . 2. 11)  is 
identical  with  the  equation  satisfied  by  amplitudes  of  the  geometrical  optics 
waves  in  the  asymptotic  theory  of  diffraction 

3.2.3  Behavior  of  the  Fields  in  the  Vicinity  of  the  Individual  Wavefronts. 
The  field  behavior  in  the  vicinity  of  the  individual  wavefronts  can  be 
obtained  by  considering  the  asymptotic  behavior  of  the  transform  of  the  field 
vectors  for  large  values  of  the  complex  frequency  s.  This  implies  a 
Tauberian  theorem  which  is  obtained  from  an  Abelian  theorem  given  by 
Bremmer  and  van  der  Pol  (Brown,  1962)  as: 

If  a  one-sided  original  jliere  time  function  f(t)  is  represented  asymp¬ 
totically  as  t  — >0+by  some  power  series  of  not  necessarily 
integral  exponents  exceeding  -1,  their  the  series  in  s,  obtained 
by  transposing  the  original  term  by  term  represents  the  image 
Qiere  Laplace  transform  F(s)  =  f(t)  ^asymptotically  as  s  -*  oo. 
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However,  here  the  inverse  of  the  theorem  must  hold  as  well  which  will 

enable  to  determine  the  asymptotic  behavior  of  the  original  time  function 
+ 

as  t  — >  O  from  the  asymptotic  representation  of  the  image  as  8  -md. 

In  order  to  justify  the  use  of  the  "Tauberian  Theorem"  which  is  the  inverse 
to  the  above  Abelian  theorem,  the  nature  of  the  time  functions  resulting 
from  the  inversion  of  the  asymptotic  representation  of  the  image  must  be 
considered. 

A  time  function  of  the  form: 


^(r,  t)  =  a  6 


t  -  t 


oo 

+  U  (t  -  t  )  ^  a  t 

o  < — v  n 

n=0 


U  —  O' 


(3.2.13) 


has  a  Laplace  transform  if  a  >  -1,  where  the  inverse  is  given  as 


-s  t 


F(r,  s)  =  e 


oo 

a  +  V' 


P(u  -  a  +  1) 


(3.2.14) 


a  >  -1 


consequently,  since  the  transform  of  the  field  vectors  is  obtained  as  a 
solution  of  a  vector  wave  equation  and  if  for  such,  in  general,  a  Laplace 
transform  exists,  it  can  be  said  that  the  transform  of  a  vector  field  must 
be  representable  in  the  form  of  Eq.  (3 . 2. 13)  in  the  vicinity  of  a  wavefront 
(for  a  rectangular  pulse  similar  holds).  With  the  above  stated  theorem, 
allowing  for  an  asymptotic  expansion  for  large  s,  it  is  reasonable  to  assert 
that  the  asymptotic  representation  of  the  inverse  of  these  transforms  can 
be  obtained  by  a  term  by  term  inversion  of  their  asymptotic  representation 
for  large  s. 

Using  the  expansion  of  Eq.(3.2.9),  it  can  be  shown  that  for  large  values 
of  s  this  equation  can  be  rewritten  as: 
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ltm 

s»l 


F(3)  ’  E  a‘>>lUm(i'  s)  exp  t8  +  2  «  -  8  +  Tm(^]  <3' 

m=u 


The  time  function  which  corresponds  to  the  mth  term  in  (3. 2. 15)  will  be  zero 

for  t  <  T°(r).  Since  T°  (r)  was  chosen  such  that  it  satisfies  the  Eikonal 
m  —  m  — 

equation,  the  equation  t  =  T  (r)  describes  the  arrival  of  the  mth  wavefront 

at  r.  The  amplitude  of  the  disturbance  associated  with  the  wavefront  is 

given  by  lim  U  (r,  s).  The  function  U  (r,  s)  satisfies  for  s  — woo: 
m  —  m  — 

s  -*oo 


2  VU  (r,  s)  •  VT°  (r)  +  U  (r,  s)  V2  T°  (r)  =  0 
m~  m  —  m—  m  — 


(3.2.16) 


which  is  the  equation  satisfied  by  the  amplitude  of  waves  determined  by 
geometrical  optics  Thus  it  is  apparent  that  in  pulse  diffraction 

problems  the  amplitudes  of  the  waves  at  the  wavefronts  are  determined  by 
geometrical  optics. 

3.2.4  Behavior  of  the  Fields  when  [t  -  T°  (r)J  >>  1. 

- - - - w- - m  —  - 

Neglecting  the  losses  (cr  — >  0),  the  inverse  transform  for  the  mth  term 
in  the  wavefront expansion(3. 2. 9) can  be  formulated  as: 


Jl  1  [f  (r. 


/•B  +i  oo 


s) 


1 

2jr  i 


U  (r,  s)  e 
m  - 


t  -  T  (r) 
m  — 


ds  (3.2.17) 


B-ico 


where  a  typical  amplitude  function  U  (r,  s),  in  general,  will  have  both  poles 

m  — 


t  -  T° 
m 


(r)J 


>>  1 


and  branch-points  in  the  region  R(s)  <  B.  The  fact  that 
necessitates  the  deformation  of  the  contour  of  integration  in (3. 2. 17)  into  the 
left  hand  plane,  where  Re  (s)  <0.  If  t  -  T^(r)  R  >>  1,  the  most  signif¬ 
icant  contribution  to  (3. 2. 17)  results  from  the  singularity  which  has  the  maximum 
real  part.  In  Fig.  3-3,  a  distribution  of  poles  and  branch  cuts  is  demonstrated 
for  the  case  of  plane  wave  diffraction  by  a  finite  convex  body.  In  this  case 


2.15) 
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FIG.  3-3:  INTEGRATION  IN  THE  COMPLEX  S  PLANE. 


C  ^  *  costf^+sin©^  x 
C  *  cos  © _  2  + sin©  x 

o  £*  & 


(3.3.7) 


x 


FIG. 3-4:  GEOMETRICAL  CONFIGURATION. 
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the  pole  along  the  imaginary  axis  at  s  =  -  iuQ  and  the  branchpoint  8  =  c^q  =  0 
are  the  significant  singularities,  where  the  contribution  from  s  =  -  iw^  yields 
the  CW  behavior  and  the  integration  around  the  branchpoint  c^  =  0  will  describe 


1 


the  deviation  of£  Im  (r,  s)  from  its  steady  state  behavior.  Depending 


upon  the  relative  value  of 


1,  the  contributions  of  the  other 

R »  1. 


{*  -  T°m 


singularities  may  have  to  be  encountered,  not  so  for 

In  fact  for  t  -  T°  (r)  R  —  1  the  exact  solution  of  F  (r,  s)  is  required, 
L  m  -J  m  *■' 

which  in  most  of  the  cases  can  only  be  approximately  solved. 


3.2.5  The  Dispersive  Effect  Due  to  the  Presence  of  Finite  Con¬ 
ductivity  or  the  Properties  of  Layered  Media 

3.2. 5.1  Finite  Conductivity. 

Although  the  constituent  parameters  ij,  e,  cr  have  been  assumed  to  be 
individually  of  a  nondispersive  nature,  the  fact  that  cr  ^  0  introduces  a  dis¬ 
persive  effect,  since 

e1  *  €(  1+—)  (3.2.18) 

es 

Assuming  a  non-conducting  propagation  media,  the  inherent  source  of  dis¬ 
persion  then  may  be  related  to  the  effects  of  the  conducting  obstacle  only. 

Applying  the  results  of  section(3. 2. 3)  toequation  3. 2. 18),  the  behavior  of  the 
reflected  wave  solution  in  the  vicinity  of  its  wavefront  is  determined  as 

lim  =  lim  e  (1  +  ~ — )=  e  (3.2.19) 

€s 

8  —*>00  8— >00 

This  indicates  that  even  a  good  conductor  acts  like  a  dielectric  in  the 
immediate  vicinity  of  the  wavefronts,  i.e.  the  reflection  coefficient  is 
determined  by  the  dielectric  constant  only.  This  property  will  be  shown 
in  the  next  section,  where  the  transforms  of  the  individual  field  vectors 
are  derived  explicitely  for  the  scattering  from  a  lossy  slab  mounted  on  a 
perfectly  conducting  planar  surface. 
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As  long  as  oblique  incidence  is  excluded,  in  general,  it  can  be  said 

that  the  dispersive  effect  of  a  conducting  media  is  solely  caused  by  the  fact 

that  a  f  0.  The  significant  parameter  in  the  assessment  of  the  magnitude 

of  the  dispersive  effect  is  the  relaxation  time  t  3  (2jr)/  (a  /  e),  where 

c 

u  =  (2jr)/  r  3  (a  /  e)  may  be  denoted  as  the  critical  radiant  frequency  of 
c  c 

the  medium.  Ifo»  1,  i.  e.  u  >>  1,  then,  in  general,  the  dispersive 

c 

effect  can  be  neglected  since  signals  decay  in  a  time  of  a  few  relaxation 

periods  t.  It  appears  that  the  dispersive  effect  cannot  be  detected  if 

u  >>  w  where  w  is  determined  by  the  high-frequency  cut  off  of 
c  max  max 

the  detecting  device.  Thus  only  for  those  materials  for  which  u>  <  u>  the 

c  o, 

operating  frequency,  dispersive  effects  are  non-negllglble. 

3.2.  5. 2  Layered  Structure 

Yet  there  still  exists  another  source  of  dispersion,  although  v,  e,  o 
are  of  non-dlsperslve  character  and  a  may  even  be  zero.  For  Instance  in 
the  case  of  oblique  incidence  on  a  grounded  dielectric  slab,  dispersion  will 
be  encountered  due  to  the  superposition  of  waves  which  results  in  a  sep- 
eratlon  of  phase  and  group  velocities,  deviating  from  the  intrinsic  velocity 
of  the  medium. 

3.3.  Pulse  Scattering  from  a  Lossy  Dielectric  Layer  Mounted 
on  a  Perfectly  Conducting  Planar  Surface. 

To  receive  more  insight  into  the  problem  of  pulse  scattering  from 
finite  coated,  convex  scatterers,  the  planar  case  will  be  treated  first  for 
both  oblique  T  E  and  T  M  plane  wave  incidence  for  the  general  case  of  a 
lossy  dielectric  layer.  For  simplicity  a  rectangular  pulse  of  pulsewldth 
T,  carrier  frequency  u  and  amplitude  fo  (r  *  0,  t  *  0)  =  1  is  chosen, 
propagating  at  an  angle  6 ,  to  the  normal  of  the  interface  as  illustrated  in 
FIG.  3-4.  The  corresponding  field  quantity,  normal  to  the  plane  of  Incidence 
(in  the  TE  -  case  the  E  -  vector,  in  the  TM  -  case  the  H  -  vector)  is  then 
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described  as: 

f0  1  <  ?/c 

(kr,  t)  «  x  Jexp  i  (k?  -w  t)  £/c  <  t  <  £  /  c +T 

to  5/c  +  T  <  t 

where  £  *  (cos  0  z  +  sin  0  y) 

f  *  (k  r.  t)  *  f J  (k)  exp  -  iu  t 
i  o 


(3 

(3 

(3 


Fj  (k,  r)  =  exp  k  £ 
w  =  i  S,  8  =  -iu 


,L 


This  time  dependent  function  F  (k  r,  t)  in  turn  can  be  expressed  in  terms 
of  an  inverse  Laplace  transform: 


i  .  $ 

F  <ks-  art 


p+i  co 
1  etS  Fj  (s, 
Jb  -i  co 


r+i  co 

ts  -(s+iu>  ) 
e  0 


(3 

(3 


jB-i 


(s  +  iu  ) 


00 


F^s.r) 

ds 


Any  other  shape  of  finite  impulse  can  be  obtained  by  convolution  properties 
which  will  not  be  considered  here. 

The  scattered  field  in  the  direction  of  reflection  is  given  as: 


f  (k  r,  t) 


Where 


1 

OB+i  co 

A 

tsf 

jl-e  " 

x  ' 

\  e 

2tt1 

jB-ioo 

(s  + 

=  R 

V  M 

)  (it) 

cu  v  c  ' 

i  k 
c 

8  is  . 

F  — ,  r  )  ds 
c  “ 


r8  (1SL 
'  c 


£  *  (  -  cos  0,2  +  sin  0  y) 

s  i  i 


0 

(3 

(3 


.3.1) 

.3.2a) 

,3.2b) 

.3.2c) 
.3. 2d) 

(3.3.3) 

1.3.4) 

.3.5a) 

•3.5b) 
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The  total  reflection  coefficient  R^'^  ie  determined  from  the  C  W 

C  W  c 

case,  where  v  denotes  either  TE  or  TM  incidence,  and  /j  whether  of  electric, 
magnetic  or  power  type.  For  proper  expansion  of  the  total  reflection  co¬ 
efficient  the  derivation  of  the  C  W  electric,  magnetic  and  power  reflection 

and  transmission  coefficients  is  thus  required  for  the 

v  w  u  w 

1/  n  Is 

behavior  at  one  interface  where  '  (  —  )  describes  the  case  of  a 

L  W  C 

grounded  dielectric  slab  and  must  be  derived  as  well. 


3.3.  1  Derivation  of  C  W  Reflection  and  Transmission  Coefficients. 
The  constituent  parameters  are  given  as: 

2  1  2 

a)  Medium  (1):  e,  c  e  ,  p  =  p  ,  a  *  0,  v  ° -  *  c 

1  o  1  o  1  *1  ^ 


(3 


b)  Medium  (2):  €2  *  *2  ( 1  + _ —  — ) '  ^2  '  ^o’  °2  *  ° 

€2  8 


N  (s)  =  N  (  1  +  UcP  2 


2  2  2  2 

k  =  *  -  s _ N  (s), 

L  2  2 


N2  =  ^2_,  W  *  ^_2_ 

e  0  e 
o  2 

2  2  m~2  ,  * 
v2  *  c  N  (s) 


(3 


d  =  depth  of  layer, 
c)  Medium  (3):  a  ^  =  oo 

3.3. 1. 1  TE  -  Incidence  for  d  =  oo:  f  *  (r,  t)  8  y"  E*  (r,  t) . 

Applying  the  boundary  conditions  of  equation  (2.  6)  the  electric  (voltage) 
TE  reflection  coefficitne  for  the  case  of  d  »  oo  is  given  as: 


EVTE  =  -x  =  Ev  = 


E 


cos  01  -  /n2  -  Sin2  01 
o _ _ 

^  cos  +  /n2  -  Sin2  ei 


(3 


.3.6a) 


.3.6b) 


3.8a) 
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and  the  magnetic  TE  reflection  coefficient  as. 

Hy TE  _  _EyTE  (3.3.8b) 

which  would  be  identical  with  the  electric  reflection  coefficient  for  the  case 
of  reversed  propagation,  i.e.  if  the  wave  is  traveling  from  the  denser  into  the 
less  denser  medium,  which  will  be  denoted  as 


Ey  TE  _  _  EyTE  _  HyTE  _  Hy  TE 


where  the  Fresnel's  law  of  reflection  is  given  as 
The  transmission  coefficients  are  given  as 


0 

sin  2 

sin  0^ 


N 


(3.3.8c) 

(3.3.9) 


2  (—  )  cos  0 

EWTE  =  Hw  TE  =  1  +  EvTE  =  —  ° 


(3.3.10) 


cos  0 


2  2 
sin  0 


HwTE=Ew  TE  =  i  +  HvTE  =  i-EvTE 


/  2  2 

j2  N  -  sin  0^ 


(~~j  cos  01  +  /n2-  sin201 


(3.3.11) 


The  power  reflection  and  transmission  coefficients  are  then  defined  as 


pv  TE  +  pwTE  =  1 


(3.3.12a) 


where 


PVTE  .  -  EVTE  ■  “vTE  and  PWTE  =  EWTE  .  HWTE 


(3.3.12b) 


3.3. 1.2  TM  -  Incidence  for  d  =  oo:  \r,  t)  =  x  H1  (r,  t). 

Applying  the  boundary  conditions  (3.2.6)  the  magnetic  (current)  TM 
reflection  coefficient  is  given  as 
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(3.3.13a) 


H 


IV  TM  =  _JL  =  kv  = 


— ”)  N2  cos  6 
M2 


l  -  ^N2-  sin201 


H 


jj2)  N2  cos  ^  +  /n2  -  sin201 


and  the  electric  TM  reflection  coefficient  as: 

EyTM  _  IV  TM 

which  again  is  identical  with  the  magnetic  reflection  coefficient  for  the  case 
of  reversed  propagation,  i.e.  if  the  wave  is  traveling  from  the  denser  into 
the  less  denser  medium,  which  will  be  denoted  as: 


Hy  TM  _  HyTM  _  EyTM  _  Ey 


TM 


The  magnetic  TM  transmission  coefficient  is  defined  as 


H,„TM  .  +  Hy  TM  E„.  TM 


2©  n2>:os«i 


W  =  1 


=  W 


EW™  =  HW  ™  =  1  +  Ev™ 


N2  cos  ^  +  /n2-  sin201 
2  ^N2-  sin201 


—2’)  N2  cos  d  +  /n2  -  sin2# 

V  1  1  1 


The  TM  power  reflection  and  transmission  coefficients  are  defined  as: 


P  TM  P  TM 

Jry'li*i+  W  =  1 


where 


(3.3.13b) 


(3.3.13c) 


(3.3.14a) 


(3.3.14b) 


(3.3.15a) 


Py  TM  HyTM  EyTM 


(3.3.15b) 
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and 

P  TM  H  TM  E  TM  (3.3.15c) 

w  =  w  •  w 

3.3. 1.  3  Derivation  of  the  Total  Reflection  Coefficient  R1^  'f*  (r,  s) 

tot 

V  u 

The  total  reflection  coefficient  R  .  (r,  s)  will  be  determined  by  two 

tot 

different  methods,  namely 

a)  The  summation  of  all  phase-delayed  partial  reflection  coefficients,  and 

b)  The  application  of  the  equivalent  transmission  line  concept. 


3.3. 1.3.1  Summation  of  Phase  Delayed  Partial  Reflection  Coefficients. 
The  reflection  coefficients  at  the  perfectly  conducting  screen  are  given 


as 


V  TE  lim  V  TE 

U  —  V 

N-m» 


H  TE  V  TE 
U  =  -  U 


(3.3.16a) 


H  TM  lim 

N-*oo 


TM 


=  +  1; 


The  delay  factor  A  is  defined  by 


A  »  exp  (  ^  ^2  ^  =  exp 


£ 


VU™ 


N(s) 


HU™ 


(3.3.16b) 


(3.3.17a) 


where 


jj*  =  d  cos  0  +  d  sin  0  tg 

2  2  2  2 


(3.3.17b) 


cos  6. 


this  however  is  exactly  the  desired  result,  namely  that  the  delay  factor 
is  related  to  the  group  velocity,  where 
„  1 


V 


pl2  j  2  cos  ©2  N  (s) 


(3.3.18a) 


'V  \  »  C  cos  6n 

gr.  )  „  2 


2/2 


N  (s) 


no  return  in  z-direction  for'' 
grazing  incidence. 


(3.3. 18b) 
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The  group  delay  then  can  be  defined  as: 

o  N  (s)  d  _  N  d 

8  "  C  cos  C  cos 


r~  w 

J  i  +  c 


g 


CJ 

1  +  c 


A  =  exp  <  -  s  T 


^-s  jl  +  ^c  |  =  exp  |-sT°(s)| 


(3.3.19a) 


(3.3.19b) 


where 


T  (s)  =  2  T 

n  g 


JTTaT 


(3.3. 19c) 


n  =  0,  1, 


A  =  exp 
n 


f8  Tn (8)  | 


(3.3. 19d) 


The  algorism  for  the  partial  pulse  returns  can  be  written  down  immediately 
from  inspection  of  figure  2: 


F8  (s) 
o 

F8  (s) 
F8  (8) 
F8  (s)  = 


=  V 


=  W-  A-U  -A  W 


=  W-A-U-A-V  -A-U'A-W 


W- A-U  -A-V  -A-U  -A-V  •  A-U  -A  -W 


=  V 


=  W-  W_-A  -U  (1) 

2  2 

=  W- W_‘A  -U  (A  U  V) 

2  2  2 
=  W-W  A  •  U  (A  U  V) 


F  (s)  = 
n 


oo 


R  (s)  =  F  (s)  + 

tot  o  a__j  n 

n  =  1 


f 


1 


W-W.-A  U 


s . .  2  r  2 

Fn  (s)  =  V  +  WW_A  U  1  +  (A  U  V J 


2  n  —  1 

(A  UV") 


.  2  ,2  I  W  W~  i 

(AUVJ  +...J  =V+  ITT? 


-  2 
A  U 


U  V  J 


V  +  A  U 


1  +  A  UV 


(3.3.20a) 
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Replacing  V,  A  and  U  by  the  proper  expressions  associated  with  the 
TE  or  TM  cases,  will  yield  the  corresponding  total  reflection  coefficient, 
where  the  reflection  coefficients  will  be  rewritten  into: 


VI.)  ^  with  „  .  V„TE  . 

1  +  A  (s; 


Rtot  «  * 


V  +  A  U  _  (  l-A(s)  -  (1+A(s)  exp  (-2s  T  (s)  ) 

2  o 

1  +  A  UV  (  1+A(s)  -  (l-A(s)  exp  (-2s  T  (s)  ) 


Lxp  (sT°(s))  -  exp  (-sT°(s)  ) 

-  A(sJ 

exp (sT°(s))  +  exp  (  -sT°(s)  ) 

|exp  (sT°(s))  -  exp  (-sT°(s)) 

+  A(s) 

exp  (sT  (s))  +  exp  (  -sT°(s) ) 

thus 


R  (8)  =  tan  hi  (s  T  (s)  )  -  A  (s) 
tan  hi  (s  T°  (s)  )  +  A  (s) 


(3. 


a)  Normal  Incidence:  A  (s)  =  N  (s) 
V 


J  %  V  TE  ,  „  4  tanh  s  T  (s) 

R  (e  •  0,  S  )  *  R  0  =  o,  s)  - - rTT^ 

tot  1  tot  1  .  ,  O.  , 


tanh 

-  N  (Sl 

tanh 

s  T°(s) 

+  N  (s) 

i  tan 

JJ  T°  (w) 

-  N  (w) 

i  tan 

u)  T°  (u>) 

+  N  (u) 

V  TE,  cos  9 

b)  TE  -  Incidence  =  A  (s)  =  _ 2_ 


_o_  \  N  (s) 


cos  0 


1 


M, 


(3. 


V  TE  ,  , 
Rtot  (s) 


cos  0^  tanhi 

a  T°  (s) 

-  N  (s)  cos  0 „ 

Ld 

cos  0  tanhi 

s  T°  (s) 

+  N (s)  cos  0 

Cd 

A* 


>  ,  VaTM.  .  /Mo\ 

c)  TM  -  Incidence  =  A  (s)  =  J 


cos  0 ^ 
cos  0, 


N  (s) 


(3 


3.20b) 


3.21a) 


.3.21b) 
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(3.3.21c) 


V  TM  .  . 
Rtot  (8) 


/cos  e2 

tanhi 

|s  T°(s) 

-  cos  Q ^ 

N(s) 

\  COB  *2 

tanhi 

|s  T°  (s) 

+  cos  Q 

N(s) 

3.3. 1.3.2  Application  of  the  Equivalent  Transmission  Line  Concept. 

The  equivalent  transmission  line  concept  is  derived  from  proper 
definition  of  equivalent  impedances  which  describe  the  properties  of  the 
structure  as  illustrated  in  figure  4.  The  equivalent  transmission  line 
input  Impedance  in  this  case  is  given  as: 

Z. 


Z.  (s) 
in 


lim 

Z3^° 


Z2  +  tanhi  I 

Lx 

0  (J 

1 

1  +  Z3  tanhi  | 

[ST° 

is] 

=  Z^  tanhi  1 8  T  (s) 


[s  T°  (s)] 


where  the  "voltage"  (here  electric)  reflection-coefficient  is  defined  as: 


vR  w  . 

Rtot  Z.  +  Z 


in 


1 


tanhi 

r  w 

CQ 

O 

H 

1 _ 

I-* 

tanhi  | 

\b  T°  (8)] 

hi 

thus  to  verify  the  result  of  the  preceding  section  it  is  sufficient  to  show 
V  v 

that  the  A  (s)  Z^  (s)  and  both  derivations  are  identical. 


a)  TE  -  case. 

Z 


1 


1 


cos  6 


1 


the  characteristic  impedance  of  medium 

(1) 


'1  = 


(3.3.22a) 


(3.3.22b) 


(3.3.23a) 
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N(s)cos  6, 


the  characteristic  impedance  of  medium 

(2)  „  . 


(3.3.23b) 


V.TE  i 
A  (s)  = 


'cos  6 


cos  6 


N  (s) 


(3.3.23c) 


which  is  identical  with  the  result  obtained  before. 


b)  TM-case 

Z  =  Z  cos  6 
1  °i  1 


where 


(3.3.24a) 


,  z  cose 
°2  N (s) 


where 


(3.3.24b) 


VA™  (s)  «(  ~- 


p  \  cos  0 
o  \  1 


2  /  TM  \  M2 


cos 


N  (s) 


(3.3.24c) 


which  also  checks  with  the  result  obtained  before.  The  relationship  be¬ 
tween  the  total  reflection  and  total  transmission  coefficients  is  identical 
to  that  derived  for  the  partial  coefficients.  Thus  the  total  reflection  coeffi¬ 
cients  always  can  be  expanded  as 


R.  f  (8)  =  F®  (s)  + 
tot  o 


F8  (s)  =  V  +  (  W-  W  -  A2  U) 


=  1  +  (A  U  V_ )  +  (  A  U  V_) . 

2 

where  A  takes  care  of  the  proper  time  delay  of  the  pulse  returns  which  is 
strictly  related  to  the  group  velocity  as  shown  by  Eqs.(3.3. 18b)  and 
(3.3.19b). 


(3. 3. 24d) 
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3.3.2  Evaluation  of  the  Inverse  Laplace  Transform. 

Having  expanded  the  C  W  reflection  coefficient  R^^s)  into  the 
proper  algorism,  the  backscattered  field  can  be  derived  by  evaluating  the 
Inverse  Laplace  Transform  of  equation  (3.3.4): 


f  S(r,  t) 


-  _i_ 
2nl 


rtB+  ioo 
J  B-ico 


e+ts 

(s  +  iu  ) 
o 


ik? 

e 


ds  . 


For  proper  evaluation,  it  is  convenient  to  introduce  anew  special  time 

variable  t  by  choosing  an  observation  point  at  (0,  0,  z,  ct  ) 
o  o 

t  =  -+ 1  .  (3.3.25) 

c  o 

Introducing  this  variable  into  (3.3.4)  and  rewriting  the  Inverse  Transform 
into  two  parts,  it  can  be  shown  that 


fS(r,  t)  =  f8  (r  6  t  )  +  f8(r,  t  ) 
—  —  l  —  o  i  o 


(3.3.26a) 


where 


(r.t  ) 
—  o 


_JL 

27T  i 


flB+ico 
J  B-ico 


exp  ^2/ o+tQ]  s 

(s  +  iu  ) 
o 


ds 


(3.3.26b) 


a  -iu0  T 

fl(r.t)  -  J-S  -°- 
—  2  —  o  2  7T  X 


joB  +  i  00 

V  B-ioo 


exp  [z/c-Ht  -T)]s 


s  +  iw 


R. 


tot 


(t) 


ds 


-iD  T 

=  -  e 


T) 


thus  for: 

t<2/c  or  t<0  : 

'  o 

Z/c  <  t  <  Z/c  +  T  or  0  <  t  <  T: 
/  _  _  /  —  o— 


(3.3.26c) 


fS(r,  t  )  =  0 
“  —  o 

f  8  (r,  t  )  =  f  8  (r,  t) 
-  —  o  —  i 
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(3.3.27) 


t  >  Z/  c  +  T  or  t  >  T:  l8  (r,  t  )  *  f  8  (r.  t  ) 

o  o  o 


t*  i  .V  tBt  -iwT  ,8.  ^ 

+  f  <£/  t  )  *  f  1  (r,  t  )  -  e  o  f,(r,  t  -T) 
2  o  1  o  l  o 


The  last  expression  rewritten  into  Fourier  Transform  terminology  becomes: 


f  (r,  t  )  ■  f 8  (r,  t  )  »  f8(r,t) 

— '  O  1^0  1^0 


1-U  (t  -T) 
o 


<: 


where  U  (t)  is  the  Heaviside  unit  step  function,  and  thus  for  t  >  T  the 
initial  response  is  cancelled. 

3.3.  2.1  Normal  Incidence 

At  normal  incidence  the  TE  and  TM  cases  cannot  be  distinguished, 
and  the  electric  reflectic  coefficient  is  given  as: 


Hv  .  -Ev 


1  -  N(s) 

1  +  N(s) 


0 


and  the  delay  factor  as: 

2 

A  *  exp 
Yielding 

RVo%  *1  ■  e2  -° 


^  o 

-2  s  T 

1  +  c 

\ 

8 

o  Nd 
T  =  — 
c 


1  -  N(s)  _  4  N(s) 

1  +  N(s)  (1+N(s) 


00 


exp 


-2  s  T 


to 

1+  c 


1 


s  J  n=l  *- 


(j~)exp(-2sToJ  1+Ic_ 


(3 


n-1  ( 


let  o  o 

T  ■  n  T 
n 


1  +  u 


c 

8 


(: 


.3.28) 


.3.29a) 


.3.29b) 


.3.29c) 


.3.29d) 


39 


then 


(  oo  ,  .  N 

12..  e  .  oV  >,  AiM 

c  l  J  1+N(b)  1+N(»)  ./^l  l  1+N(«) 


n-1 


exp 


-2  s  T 


n 


3.3. 2. 1.1  Lossless  Case:  N  (s)  °  N.  T°  (a)  c  T° 

In  this  case  the  Inverse  Laplace  transform  can  be  evaluated 
straight  forwardly,  since  only  the  pole  at  s  s  -  i  exists,  and 


i8  V  T  Tn  =  y 


Rn  U  (t  -  2n  T  ) 
o 


(3.3.30) 


1-U  (t  -2n  T  T) 
o  ~ 


E  (3.3.31a) 
o 


where 

E 

—  o 


=  ^  1  dxp  (-  -  ? 1  ) 
c 


R 

o 


1  -  N 

1  +  N 


R  = 
n 


4  N 

(1+N)* 


(3.3.31b) 


(3.3.31c) 


o 

Thus  if  2  T  >  T,  the  pulse  returns  are  separated,  in  Fig.  3-5  the  cases 
N  *  2,  2  T°  =  2T  and  T  =  2  0  T°  are  plotted. 


A  meaningful  application  of  pulse  scattering  to  the  problem  of  inverse 

scattering  requires  an  operation  with  pulses  of  extremely  small  tlme- 

o 

width,  or  of  exceedingly  large  timewidth  where  T  <  2  T  .  Since  it  must 

be  assumed  that  T°  is  extremely  small  and  the  information  from  which  the 

o 

electric  properties  can  be  determined  are  the  time  delay  factor  T  of  the 

layer  as  well  as  the  amplitude  ratios  of  the  partial  pulse  returns,  these 

quantities  must  be  distinctly  measureable.  The  refractive  index  of  the  layer 

in  the  lossless  case  can  be  determined  as  the  average  value  over  the  second 
th 

to  the  m  time  delayed  pulse  returns,  excluding  the  first  return  as 
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FIG.  3-5:  PULSE  SCATTERING  FROM  A  GROUNDED  LOSSLESS 
DIELECTRIC  SLAB  (  N  =  2,  T°  =  T). 
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N  = 


(3.3.32a) 


1-  en  +  1,  n 

1-1  e  +1,  n 
J  n 


n=l 


in 


c  +  1,  n 
n 


R  +  1 
n 

R 


1  -  N 

1  +  N 


(: 


and  the  depth  of  the  layer  then  can  be  determined  from  the  time  delay 


factor  T  as 


o 

T  C 

N 


0 


CJ  I/O  W  / 

3.3.  2.  1.2  Lossy  Case:  N(s)  c  N(l  +“)  ,  T°  (s)  c  T°  (l+~)  ' 


In  this  case  in  addition  to  the  pole  at  s  *  i  u  a  branch  cut  will  be 

o 

introduced  which  degenerates  into  a  pole  at  s  *  -  ^  f°r  slightly 

lossy  case,  where  the  following  approximation  may  be  implied: 

2 


N(s)  =  N  ( 1  +  *^c  ^  l!2  - 


N  1  + 


2s  J 


l 

8 


w 


In  general,  the  partial  reflection  coefficients  for  the  lossy  case  at  normal 
incidence  are  given  as 


Ev  _  1  -  N(s) 

1  +  N(s) 


2  N 

1  +  N 


1-N 
1  +  N 


c 


(3.3 


.3.32b) 


.3.32c) 


.3.33) 


.34a) 
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1-N  ,2  N 

R  c  -  +  - 

o  1+N  1+N 


s 


(3.3.34b) 


(3.3.34c) 


(3.3.35a) 


(3.3.35b) 
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(3.3.35c) 


bnYh“c^2 


U 

+  N  1  +  c 


1/2 


^  exp ' 


-48  T 


R 


exp 


-2n  s  T 


By  inspection  of  (3.  35)  it  follows  that  in  the  vicinity  of  the  wavefronts 

of  the  partial  pulse  returns,  i.  e.  for  (t  -  2ir  T  )~0,  the  behaviour  is  iden- 

o  o 

tical  to  that  of  the  lossless  case  according  to  the  Tauberian  theorem  of  section 

2.  3.  In  general,  the  exact  evaluation  is  tedious  ,  since  two  branch  points 

are  introduced  at  s,  ■  0  and  s,  =  -u  which  for  a  small  loss  approx- 

1  2  °'  l 

imation  degenerate  into  a  pole  at  S  *  -  u  ,  since  N  (s)  — N  (l  + _ c_  “), 

P  Si  g  o  8 

N  +  1 

In  this  case  the  integrals  can  be  evaluated  straight  forwardly  and  will  be 
treated  next: 

Small  Losses  :  N  (s)  N  (  1  +  ^c  1^  ) 

2  s 

Applying  the  approximation  of  eq.  (3.  33),  the  partial  reflection 
coefficients  can  be  rewritten  into: 


R 


_  1  -  N  (s) 
1  +  N  (s) 


1  -  N 

1  +  N 


s  +  a 

s  +  b 


u 


a  =  • 


Ri*-<7?&2  <”‘p-2sT°<8> 


c/2  N 

N  -  1 

U)  i 

b=_c/2N_ 
N  +  1 

c=Uc/2 


(3 


4  N  „  s  (s  +  c) 

TTTn)2  (7Tb p  exp 


-2  T°  (s  +  uc/  2) 


(3 


(3.3.35d) 


.3.36a) 


.3.36b) 
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R 


4  N 


n  (1  +  N) 


s  (s  +  c)  2 
(s  +  b) 


n-1 

-O 

exp  -  n  u  T 
c 

o 

exp  -  2n  T  8 


(3.3.36c) 


The  total  reflection  coefficient  is  thus 


R  (s,  u  ,  T  ,  N  )  = 
tot  c  o 


1-N 

1+N 


given  as 

s+a  -  4  N 
s+b  (1+N)  2 


8(8+0) 

fe+b) 2 


oo  / ,  \  n-1  /*  \ 

E(^~)  fe-) 


n-1 


exp  (  -n  u  T°) 
c 


dxp(-2nsT  )  (3.3.37) 


n=l 


The  integration  in  the  complex  s-plane  thus  results  simply  in  the 

evaluation  of  the  pole  at  s  *  -  1  cj  and  about  the  n^  order  pole  at 

o 

(J 

s  =  -  c  N  as  illustrated  in  Fig.  3-6,  there  the  inverse  transform 

. P  2  N+l 
is 

is  given  by: 


A 


s  » b 


8=-iu  X 
o 


FIG.  3-6:  DEFORMATION  OF  CONTOUR  FOR  THE  SLIGHTLY 

LOSSY  CASE  N(s)  —  N  (1  +  u  1  ). 

c  ~ 


is  given  by: 
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R  (s,  w  ,  T  ,  N) 
tot  CO 


(3.3.38) 


*1  {r->  'o' 


3+iCD 

°B-1  oo 


s(t  +2/c) 
e  o 


(s  +  i  u  ) 
o 


The  spatial  time  variable  t  then  is  given  by  t  *t-2/c  +  2nT°  for  the 
th  o  o 

n  partial  pulse  return.  The  inverse  transforms  in  this  case  yield: 

(See  Table) 

Immediately  it  is  recognized  that  at  the  trailing  end  of  the  pulse 
(T  >>  2T°),  the  transient  contributions  only  will  be  present  if  b  is 
small.  The  complete  solution  can  be  given  as: 


f8  (t  ,  2  n  T°)  »  R  (A  e-iWto  +  B  e"*5  *o)  U  (t  ) 


o  o 


1-U  (t  -T) 
o 


00 


R 


n=i 


n 


-iwt  ,  2  n  -bt 

A  e  o  +  (B  +  B  t  +  B  t  + . B  t)e  o 

n  n  n,  n.  nn 

_  o  1  2 


U  (t  -  2n  T  ) 
o 


1-U(t  -  T  -  2n  T  ) 
o 


(a 


where 


R 


1-N 

1+N 


4  N  (1  -  N) 


n-1 


n 


(1  +  N) 


n+1 


exp  (-n  u  T  ) 
c 


(a 

(a 


and  the  A  ,  B  may  be  determined  by  Heaviside's  expansion  theorem, 
n  nn 


For  N  >>  1,  a  =  b  =  c  «  ^c  ,  it  can  be  seen  that  3.3.37  simplifies  to 

2  ^ 
t,  /  ,v  1-N  4  N  s  \  1-N  n 

Rtot  S'WC,  C  1+N  ”  (l+N)2  S-HJ,  1+N 


c/ 2  n=l 


exp 


-2n  T 


Vt  +ei 


(a 


.3.39d) 

. 3. 39e) 

.3. 39f) 


.3.40) 
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in  which  case  the  inverse  transform  of  (3.3.39a)  can  be  applied,  yielding  with 
(3. 3. 39e,  f): 

-  iut 


fS  (  t  ,  2n  T°,  N»  1)  =  R  e 
o  o 


u(t  ) 

o 


1-U  (t  -  T) 
o 


] 


-iut 


U)  CD 

"  C  t  H 


+  A  e 
o 


+  B  e  2 
o 


°  y,  R  U  (t  -2n T°)  (T-U(t  -T-2nT°) 

J  n^r  n  0  L  0  J 


Before  the  question  will  be  answered,  whether  the  parameters  N,  p,  d  of  the 
layer  can  be  uniquely  determined  from  pulse  scattering  methods,  the  case 
of  large  losses  will  be  treated  briefly. 


Large  Losses:  N  (s)  =  N  > 


u 


1  + 


c/2  1/3 


The  partial  reflection  coefficients  for  large  losses  are  defined  by 
Eq.  (3. 3. 35)  where  it  was  shown  that  the  behaviour  in  the  vicinity  is  aside  from 

/  0  V 

the  material  attenuation  factor  exp  (-n  T  u>  )  identical  to  that  of  the  lossless 

c 

case.  An  exact  evaluation  of  the  I.  L.  T.  in  general  will  imply  convolution 

methods,  where  the  pole  at  s  =  i  u>  will  be  extracted  and  the  I.  L.  T.  then  must 

c 

be  solved  for  the  case  of  delta  pulse  incidence  and  the  following  types  of 
integrals  must  be  evaluated  in  the  complex  s-plane: 

oo 


/  (r,  t  )  - 
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i_r 

2*i  I 

tfB-i 


St 

^C_  1/2 

e 

J-N(l+  s  )  

ds 

oo 


1  +  N  (  1  +^c_) 
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n+1 


(3. 


and 


A  A  XT 

rs  /  .  ,  v4N 

f  (£.  t  )  *  ~ir~. 
-n  o  2;ri 


,B+i  oo 


St 

1  N(l,Wc  >1/2 

n  -  1 

e 

1  - N  \l  +  ; 

 s  J 

l+N  (  1  +  ^c  ) 

_  s 

n+1 

a  A ,  V* 
8 


exp 


2n  T  ( 
o 


ds 


(3 


.3.41) 


3.42a) 


.3.  42b) 
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Introducing  the  following  transformation  U  =  u*  it  can  be  seen  that  both 

c 

types  of  integrals  have  branch  points  at  U  =  0  and  U  *  -1,  where  3.3.47a,  b 
become:  , 


f8  (r,  p  >  u  t)  ■ 


A 

y  u 


.B  +i  oo 


2  »i 


.  -  2n  T  u  u  (u+1)^2 
*  e  o  c  du 


0 

eUp 

riT- n  (u + 1) 

du  (3.3. 

0 

[uV  N  (  u  +  1)  ^2 

-e'+i. 

0  7 

f 

eUP  (u+l)V2  UV2  | 

J  X-i  00  ^U^2  +  N  (u+1)  ^2J  n+1 

n-1 


(3.3.43b) 


The  behaviour  in  the  complex  u  -  plane  is  illustrated  in  Fig.  3-7: 


FIG.  3-7:  DEFORMATION  OF  CONTOUR. 

Introducing  the  new  variable  U  *  +  1  ir,  where  the  positive  sign 

corresponds  to  the  top  sheet,  it  can  be  shown  that  by  adding  the  contributions 
of  both  sides  around  the  branchcut,  3.3.43a  becomes: 
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A  -8 

y  -f  (r,  P-u, 


.)  .  VL  f 

211  Jo 


-4  i  v^2  (1  -  v)^2  e VP  dv 


v  +  ( 1  -  v)  N 
-vp 


2  u  P  1/2/,  v  1/2  -vp 

*  - _ c_  |  (1  -  v)  _  e  dv 


which  for  N  *  1  only  is  tabulated,  then 


u 


c  t 


A  *  fv  8  1  **  2 

x  .  f  (r,  p=u  t)  =  -  -  7*  e 

o  c  t  t 


u 


h  ("2~  '  u(t) 


differentiation  with  respect  to  b  =  2  T  w  : 

c 


(3.3.44a) 


(3.3.44b) 


Also  the  contribution  to  the  second  pulse  return  can  be  evaluated  for  N=1  , 
using  the  transform  pair  5.  6.  45  ofjBateman  J  by  proper  repeated 


a  ,  s  ,  .  16 

y ’In  (r,  P=u  t)=  -  — 

1 1  C  /  n  _Oi  o 

u  (t  +  2  T  ) 
c 


(t  -  2  T°)  e~l/2V  I2(“)t  •  U(t-2T°)  (3.3.45a) 


Similarly  the  consequent  pulse  returns  can  be  evaluated  using  convolution 
properties  and  the  transform  pair  (  5.  6. 22)  ofjjJateman  J  The  exact 
results  will  not  be  given  here,  but  it  can  be  concluded  that  for 
t>  2n  T°: 

-  1/2  w  t 

$‘iB  (r,  p-w  C  1  (wc  t)  U(t-2nT°)  (3.3.45b) 

ni  c  2 

g 

In  Fig.  3-8the  results  of  x  f  (r,  p«w  t)  are  plotted  for  Nal 

«1  c 

It  can  be  immediately  recognized  that  the  higher  order  pulse  returns  decay 
exceedingly  fast. 
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FIG.  3-8:  x  -f8  ,  (r,p  =  ut) 
—  nl  —1  c 


1_ 

exp  (  -  2  u  t ) 

_ c 

n  +  1 
t 


n+1 
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For  N  ^  1  no  closed  form  solutions  exist  for  any  of  3. 3. 43.  Yet  enough  information 
has  been  obtained  to  derive  conclusions  about  the  possibility  of  identifying 
the  parameters  N,  u  ,  T,  from  which  e  ,  a  _  and  d  may  be  determined. 

C  £.  C 

If  the  wavefront  amplitudes  are  distinctly  determinable  and  the 

o 

pulse  returns  separated,  i.  e.  T<2T  ,  then  the  ratio 


o 

can  be  determined  and  T  from  the  seperation  of  the  individual  pulse 

returns.  The  attenuation  factor  exp  (-u  T°)  must  be  determined  from  the 

c 

oscillogram  of  the  individual  partial  pulse  returns.  Thus  for  lossy  coating 
the  pulse  scattering  method  seems  not  to  be  the  best  approach.  In  figure 
the  pulse  returns  for  N  «  const.,  d  =  constare  plotted  for  various  values 
of  a  to  illustrate  this  point. 


3.3. 2. 2  Oblique  Incidence. 

At  oblique  incidence  the  TE  and  the  TM -cases  must  be  treated  se¬ 
parately,  however  as  long  as  o  ^  =0  and  N(s)  =  N,  the  evaluation  of  the  inverse 
transforms  is  identical  to  that  of  normal  incidence  on  a  lossless  structure 
For  the  TE-case  the  reflection  coefficient  are  given  by  Eqs.(3. 3. 8)to(3. 3. 11) 
and  (3.3.2  lb)  and  the  delay  time  of  the  partial  pulse  returns  by  (3.3. 1 9a)  which 
is  the  same  for  the  TM-case.  For  the  TM-case  the  reflection  and  trans¬ 
mission  coefficients  are  given  by  Eqs.(3.3. 13)to(3.3. 15)  and(3. 3. 21c). 

Since  the  plane  wave  concept  was  applied,  it  is  arbitrary  where 

the  observation  point  is  located,  thus  let  R  (0, 0,  Z  )  and  r  «  Z  cos  6 , 

o  o 


R  U  (t  -  2n  T°) 
n  o 

TE,  TM 


see  Fig.  3-4.  The  solution  then  is  given  as: 

oo 


f8  (r,  t  ,  6  ) 
—  o  1 


TE,  TM  n 
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[l-U<t  -  2n  T°  _t"J 


.46) 
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V 


(r,  t  ) 
o 


FIG.  3-9:  PULSE  SCATTERING  FROM  A  LOSSY  GROUNDED  DIELECTRIC 

SLAB  (  N  !  2,  T°  =  T,  u  ^  0). 

c 
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One  special  case  is  of  interest,  namely  the  TM-case  at  critical  incidence, 

i.  e.  if  tan  6  *  N,  (6  —  Brewster  angle)  which  will  be  discussed  in 

Br  Br 

more  detail. 

3.3.2.3. Behaviour  at  the  Brewster  Angle  for  TM-Incldence . 

T  H  TM  E  TM  n  ^  E  TM  ,  H...TM.  111118 

In  this  case  V  *  -  V  *  0 ,  and  W  *  1  w 

it  follows  from  (3. 20)  that  there  only  will  exist  one  partial  pulse  return, 
namely  the  second  pulse  return  or  the  first  internally  reflected  one. 


where 


„  H,_TM  . 
R  1  V  *  0 
o 


V«W™  U  VM,A2.-^(-12T°U) 


R  (n>  1)  =  0 
n 


and 


TM  .  . 

r.  .  ( e  *  b  ) 

tot  1  1_ 

Br 


i  tan  (  T°  u)  -  1  ,  .  „  _o  , 

- - — 1 -  =  -  exp  (-  i  2  T  u) 

i  tan  (  T  u)  -  1 


(3.3.48) 

(3.3.49) 


This  particular  behaviour  that  at  Brewster  angle  incidence  only  one  pulse 
return  will  be  observed,  may  be  used  to  determine  the  refractive  index 
if  the  angle  of  incidence  is  known  precisely.  However  from  this  measure¬ 
ment  alone  the  depth  of  the  layer  cannot  be  determined.  To  find  some  in¬ 
formation  about  the  depth  of  the  layer,  the  angle  of  Incidence  has  to  be  changed 
such  that  the  delay  time  T°  ( Q  ^)  can  be  determined. 

In  this  case  a  £  0,  the  behaviour  becomes  too  complex  and  pulse 
scattering  at  oblique  incidence  is  not  useful.  Only  in  the  case  of  small  losses 
and  at  Brewster  angle  incidence  satisfactory  results  for  the  determination 
of  the  refractive  index  may  be  obtained. 
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IV 


ANALYSIS  OF  THE  COMPUTATIONAL  RESULTS  FOR 
THE  APPLICATION  OF  THE  BOUNDARY  CONDITIONS 


etxet 


=  0,  AND 


E 

—  II 


-  E 
|—  s 


=  0. 


4.1  Introduction. 


In  the  following  outlines,  further  computational  results  will  be  presented 

in  addition  to  those  given  in  Quarterly  Report  No.  8579-3-Q  (Weston  et  al,  1967) 

where  the  application  of  these  two  boundary  conditions  to  the  test  case  of  a 

perfectly  conducting  sphere  of  electric  measure  ka  =  2  was  discussed.  There 

the  expressions  for  E.  (incident  field),  E  (scattered  field),  and  E  _  =  E.  +  E 

—  i  — s  — T  —  i  — s 

are  defined  as  well  as  the  properties  of  the  two  boundary  conditions. 


In  particular,  it  was  shown  that 


f|E  . 

-  |e 

LI"1 

|  —  s 

} 


possesses  only  one  minimum 


apart  from  a  limited  conical  section  about  the  focal  point  in  the  shadow  region. 
Because  of  the  vector  nature,  this  condition  proves  to  be  rather  dependent 
upon  local  polarization  properties  and  in  general  will  not  yield  the  proper 
surface  locus,  apart  from  a  limited  conical  section  about  the  specular  point. 

The  boundary  condition  E^,  x  E'^,  ;  in  contrast,  is  applicable  far  into  the 
shadow  region  with  a  very  pronounced  minimum  at  the  proper  point  which  lies 
on  the  surface  of  the  scatterer.  The  application  of  this  condition  is  however 
limited  due  to  the  fact  that  in  addition  to  the  proper  surface  locus,  a  family  of 
concentric  surface  loci  (hyperboloids  of  revolution)  will  be  obtained,  though  of 
less  pronounced  minima  (see  Figs.  3-3a  to  3-3d  of  8579-3-Q)  which  were  found 
to  lie  exterior  to  the  proper  surface  locus.  The  resulting  minima  in  specular 
direction  (6  =  180°)  are  spaced  at  equidistant  intervals  of  k  A  R  =  7r/2.  Finally 
it  was  pointed  out  that  both  boundary  conditions  have  to  be  employed  simulta¬ 
neously  to  precisely  determine  the  point  which  may  lie  on  the  surface  of  the 
scatterer. 

These  properties  of  how  to  employ  the  boundary  conditions  properly, 
required  the  computation  of  further  test  cases,  including  spherical  test  cases 
of  larger  electric  measure  as  well  as  prolate  spheroids.  The  results  obtained 
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are  presented  in  the  following. 

4.2  Spherical  Test  Cases. 

Employing  the  notation  of  Stratton  ( 1941  ),  the  expansion  coefficients  for 
spheres  of  electric  measure  ka  =  5  and  ka  =  10  were  computed  as 


;s  ,.,n  n(n+l) 

i  =  (l)  — - -  a 

oln  2  n 


,.,n  n(n+l) 

(i)  — o - 


in(P> 


hU)<p) 

n 


p  =  ka 


(4.  1) 


’eln 


,.,n+l  n(  n+ 1 )  , 

(i)  - - -  b 

2  n 


,..n+l  n(n+l) 

(i)  — o - 


<pjn<p» 

n 


p  -  ka  .  (4.2) 


In  Tables  1  and  2,  the  coefficients  for  the  particular  test  samples  are  given. 

Applying  the  two  boundary  conditions  E  ^  x  =  0,  and  ||  E  .|  -  |Eg|  ►  =  0 
to  the  near  field  expansions,  computational  results  are  obtained  (see  Figs. 

4- la  to  4-lg)  which  are  very  similar  to  those  of  the  spherical  test  case  ka  =  2 
as  presented  in  Figs.  3-3a  to  3-3d  of  8579-3-Q. 

Because  of  the  symmetry  properties,  the  obtained  results  for  azimuthal 
dependence  0,  0  +  7 r  are  identical  and  it  is  sufficient  to  plot  the  resulting 
curves  for  one  0  quadrant  over  9  -  180°  (-n  7.5°)  0°,  where  in  Figs.  4- la 
to  4-lg  the  results  for  both  ka  =  5  and  ka  =  10  are  plotted.  It  can  be  seen  that 
except  for  additional  pseudo  surface  loci  within  the  proper  surface  locus  of 
the  associated  sphere,  the  distribution  of  the  surface  loci  of  either  condition 
is  identical  in  character  for  all  three  cases  ka  =  2,  5,  10.  The  presence  of 
the  additional  pseudo  loci  within  the  proper  surface  locus  of  the  particular 
sphere  in  question  required  further  investigation.  It  was  then  found  that  the 
appearance  of  these  additional  surface  loci  depends  upon  the  number  of  expansion 
terms  which  are  used  in  the  near  field  expansion.  This  is  illustrated  in 


Figs.  4-2a  to  4-2f,  where  the  values  of 


{Stx£t} 


and  -  E 


E 

—  s 


are 
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TABLE  4-1:  EXPANSION  COEFFICIENTS  FOR  THE  SPHERICAL 
TEST  CASE  ka  =  5. 


n 

B*{W 

ta  foln} 

Ee  (Seln) 

,m  (feln} 

i 

0.4124466 

-0.2173557 

-0.7760094 

0.4169158 

2 

1.200139 

-1.469722 

1.449289 

-1.886731 

3 

-0.4013607 

5.973030 

0. 1973786E-02 

0. 1088063 

4 

-5.010763 

-4.999988 

4.409093 

2.642056 

5 

4.499645 

-1.499734 

-3.963083 

6.613639 

6 

0.1782615 

1.926581 

-3.339949 

-0.5453657 

7 

-0.4877635 

0.8499482E-02 

0. 1734405E-01 

-0.6966581 

8 

-0. 1805431E-03 

-0. 8061960E-01 

0.1020575 

0.  2893283E-03 

9 

0. 9469407E-02 

-0. 1992659E-05 

-0. 2831584E-05 

0. 1128810E-01 

10 

0. 1284376E-07 

0. 8404801E-03 

-0. 9685428E-03 

-0. 1705591E-07 

n 

TABLE  4-2:  EXPANSION  COEFFICIENTS  FOR  THE 
TEST  CASE  ka  =  10. 

Re  {foln)  Im  1‘oln)  Ee  <feln} 

SPHERICAL 

Im  (geln) 

i 

-0.4878377 

-0.6096100 

-0.3913659 

-0.4880560 

2 

1.767869 

1.475888 

-1.474155 

-1.222750 

3 

2.121697 

0.8790562 

5.150302 

2.091937 

4 

-9.997529 

0.1571675 

-0.2891902 

0. 8370103E-02 

5 

6.574667 

-3.891155 

-10.79425 

6.737797 

6 

3.204104 

-7.551152 

6.915772 

-18.40077 

7 

9.203276 

24.54987 

4.963214 

10.69282 

8 

-32.51425 

-10.64596 

5.994395 

1.027457 

9 

22.34951 

-19.90205 

-11.65034 

19.71129 

10 

6.768880 

18.06850 

-24.09158 

-14.23931 

11 

-9.214209 

1.312489 

3.855611 

-15.47917 

12 

-0.1428198 

-3.334598 

4.980965 

0.3193848 

13 

0.8995424 

-0.  8892920E-02 

-0. 1527055E-01 

1.178722 

14 

0. 3389343E-03 

0.1886478 

-0.2287811 

-0.4984860E-03 

15 

-0. 3196965E-01 

0.  8517152E-05 

0. 1147929E-04 

-0. 3711489E-01 

16 

-0. 1493902E-06 

-0. 4507445E-02 

0. 5094764E-02 

0. 1908575E-06 

17 

0. 5398010E-03 

-0.1904478E-08 

-0. 2346889E-08 

0. 5992278E-03 

18 

0. 1818520E-10 

0. 5576441E-04 

-0.6109958E-04 

-0. 2183134E-10 

19 

-0. 5030335E-05 

0. 1331804E-12 

0. 1567326E-12 

-0. 5457031E-05 

20 

-0. 7627612E-15 

-0.4002248E-06 

0. 4307819E-06 

0. 8836814E-15 

57 


0  =  0° 


FIG.  4- la:  RESULTING  SURFACE  LOCI  FOR  THE  SPHERICAL 
TEST  CASES  ka  =  5.  10.  0  =  0°. 
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0  =  15° 


- {ls,|  ■  l-s l } 

- {"I-T  x  -tI  2  } 


FIG.  4-lb:  RESULTING  SURFACE  LOCI  FOR  THE  SPHERICAL 
TEST  CASES  ka  =  5,  10.  0  =  15°. 
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ka  5 
N  =  10 


0  =  30° 


ka  =  10 


6  =  180° 


- {|E,I  -  ISj) 

-  {”  l-T  X  -tI  2} 


FIG.  4-lc:  RESULTING  SURFACE  LOCI  FOR  THE  SPHERICAL 
TEST  CASES  ka  =  5,  10.  0  =  30°. 


6  =  90° 


60 


ka  =  5 
N  =  10 


0  =  45° 


0  =  0 


ka  =  10 
N  =  20 


0  =  90 


- {|S,  I  -  I  — gl} 


-  |ETxE='  |2'- 


FIG.  4- Id:  RESULTING  SURFACE  LOCI  FOR  THE  SPHERICAL 
TEST  CASES  ka  =  5.  10.  0  =  45°. 
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0  =  60° 

ka  =  5  ka  =  10 

N  =  10  0  =  0°  N  =  20 


I  r 20 


9  =  180° 


- {|E,I  -  |E„|} 

-  {■  I E  T  x  -  T  I  2} 

FIG.  4-le:  RESULTING  SURFACE  LOCI  FOR  THE  SPHERICAL 
TEST  CASES  ka  =  5,  10.  0  =  60°. 


9  =  90° 
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ka  =  5 
N  =  10 


ka  =  10 
N  =  20 


0  =  75° 
0=0° 


e  =  iso0 


- {|E,|  -  llj} 

- {-lsT*rTi2} 

FIG.  4- If:  RESULTING  SURFACE  LOCI  FOR  THE  SPHERICAL 
TEST  CASES  ka  =  5,  10.  0  =  75°. 
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0  =  90° 

ka  -  5  ka  =  10 

N  =  10  0  =  0°  N  =  20 


FIG.  4-lg:  RESULTING  SURFACE  LOCI  FOR  THE  SPHERICAL 
TEST  CASES  ka  =  5,  10.  0  =  90°. 


9  =  90° 
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FIG. 


OF  EXPANSION  TERMS  n=l.  2,  3,  4,  5,6.  0  =  0°  9  =  135°  . 
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FIG.  4-2b:  PLOT  OF  ^(  | E  4 1  -  |Eg||}  VERSUS  THE  RADIANT  VECTOR  FOR 
THE  SPHERICAL  TEST  CASE  ka=2,  FOR  A  VARYING  NUMBER 
OF  EXPANSION  TERMS  n=l,  2,  3,  4,  5,6.  0  =  0°,  0  =  135°  . 
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OF  EXPANSION  TERMS  n=3,  5,  6,  10.  0  =  0°,  9  =  135°. 
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FIG.  4-2d:  PLOT  OF  (|jE  1 1  -  |ES  ||)  VERSUS  THE  RADIANT  VECTOR  X=kR 

FOR  THE  SPHERICAL  TEST  CASE  ka=5,  FOR  A  VARYING  NUMBER 
OF  EXPANSION  TERMS  n=  3,  5,  6, 10.  0  =  0°,  9  =  135°. 
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OF  EXPANSION  TERMS  n  =  8,  10,  11,  12,  14,  20.  0  =  0°,  9  =  135°  . 
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FIG.  4-2f:  PLOT  OF  HlEil  -  |Es||j  VERSUS  THE  RADIANT  VECTOR  X=kR 

FOR  THE  SPHERICAL  TEST  CASE  ka  =  10, FOR  A  VARYING  NUMBER 
OF  EXPANSION  TERMS  n  =  8,  10,  12,  14,20.  0  =  0°,  6  =  135°. 
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plotted  respectively  versus  the  radial  dependence  for  0  =  0,  6  =  172.5°. 
Inspecting  the  curves  closely,  it  can  be  seen  that  these  undesired  pseudo  surface 
loci  of  the  condition  will  not  appear  if  the  number  of  expansion 


terms  n  is  equal  to  the  electric  measure  ka  or  if 


n  <  [ka+(ka)*^J 


(4.3) 


whereas  with  increasing  deviation  of  n  from  this  value  such  additional  minima 
will  be  present  and  for  large  n  ,  i.e.  n  >  £ka+  (ka)1^3J  these  minima 
again  show  equidistant  spacing  of  k  AR  -  i/2  .  For  constant  dimensions  of 
the  scatterer,  the  number  of  such  minima  will  increase  with  decreasing 
operational  frequency,  where  the  innermost  is  exterior  to  the  spherical  caustic. 
Furthermore,  it  can  be  seen  that  these  additional  hyperbolic  pseudo  loci 
show  minima  which  are  relatively  small  compared  to  the  pronounced  dip  of 
the  proper  surface  locus  and  disappear  in  the  vicinity  of  the  latter,  then  appearing 
again  exterior  to  the  pronounced  dip  of  the  proper  surface  locus  as  indicated 


in  Fig.  4- la  by  the  curves  C  and  C 


Pin 


This  phenomenon  of  additional 


ext 


surface  loci  for  the  case  where  the  employed  number  of  expansion  terms  is 
much  larger  than  the  electric  measure  (4.3),  may  be  explained  by  the  fact 
that  in  the  geometrical  optics  limit  these  interior  surface  loci  are  caused  by 
rays  emanating  from  the  caustic  of  the  sphere.  This  explanation,  however,  does 
not  hold  for  the  case  where  n  is  smaller  than  ka,  these  minima  are  due  to 
zeroes  of  lower  order  Bessel  functions  of  the  vector  spherical  harmonics. 

Inspecting  the  behaviour  of  jjE  J  -  |E  J  |  for  varying  number  of  expansion 

terms,  it  can  be  seen  that  for  n<  Jka  +  (ka)1^'^  the  obtained  minima  will 
oscillate  about  the  proper  locus  and  only  for  n  >  £ka  the  minima 

will  occur  at  the  proper  surface  locus,  however,  depending  upon  local  polarization 
properties,  only  within  a  limited  cone  about  specular  direction  of  apex  angle 
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0^—  30  .  Except  for  the  conical  region  about  the  focal  point  in  the  shadow 
region  of  apex  angle  0 ^  60°  this  condition  will  only  show  one  minimum 
in  any  other  radial  direction. 

Summarizing  the  results,  the  following  can  be  concluded  for  the  simultaneous 
application  of  both  boundary  conditions  for  the  precise  determination  of  the 
proper  surface  locus: 

1)  First  apply  the  boundary  condition  ^|eJ  -  |Ej|  ,  starting 
at  the  lower  bound  k  R^  ^  .  5  ,  and  search  for  the  first  point 

Y1q  where  {l-il  -  |Ej|  reverses  signs  or  |  ||E.|  -  |E  jj  |  <  1 
with  a  searching  increment  of  approximately  k  AR  —  .25. 

2)  Apply  the  boundary  condition  ^E^,  x  E*j, |  about  this  point 


within  the  range 


Y10-W2  <  ^  <  Y10+„/2 

employing  a  rather  small  searching  increment  of  .01  <  k  1 

to  eliminate  additional  minima  and  search  for  the  proper  point  which 
shows  the  most  pronounced  minimum. 

3)  In  addition,  it  is  advisable  to  employ  the  procedure  at  two  or  more 
different  operational  frequencies  and  plot  the  points  of  successive 
minima  along  one  ray  versus  R  ,  since  then  the  loci  of  additional 
minima  will  be  shifting  whereas  the  proper  one  will  be  independent 
of  the  origin.  In  Fig.  4-3  these  properties  are  illustrated. 

4.3  Prolate  Spheroidal  Test  Cases. 

The  computation  of  the  near  field  was  again  based  upon  an  expansion  into 
vector  spherical  wave  functions  (Stratton,  1941),  where  the  approximate  expansion 
coefficients  were  provided  by  P.C.  Waterman  (  1965).  The  following  three 
test  cases  are  chosen 

a)  ka  =  1,  a/b  =  2,  n  =  6 

b)  ka  =  10,  a/b  =  5,  n  =  16  . 

c)  ka  =  20,  a/b  =2,  n  =  28  . 
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FIG.  4-3:  RESULTING  SURFACE  LOCI  FOR  A  SPHERE  OF  RADIUS 
a  ,  DETERMINED  FOR  k  =  .5,  k  =  1  AND  PLOTTED 

VERSUS  THE  RADIANT  DEPENDENCE  R  TO  ILLUSTRATE 
THE  SHIFTING  OF  ADDITIONAL  SURFACE  LOCI  (C.  , 


6  =  90° 
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The  computational  results  are  presented  in  Figures  4-4a  through  4-7d. 

Inspecting  these  figures  it  can  be  shown  that  within  a  conical  section 
about  specular  direction  of  apex  angle  0^  —  30°  ,  rather  satisfactory  results 
are  obtained.  The  deviation  from  the  actual  (proper)  locus  in  the  region  exterior 
to  this  domain  have  been  expected  by  theory.  Comparing  the  results  of  the  two 
boundary  conditions  it  can  be  said  that  within  a  narrow  conical  section  about 
specular  direction  the  condition  x  E*  j  =  0  is  superior  to  {N  -  lEsl}  =  0- 

where  it  was  found  that  the  resulting  locus  for  |  E  -  |£s|j  lies  interior 
to  both  the  proper  locus  and  that  for  |e^x  E*,  •  =  0 .  It  was  shown  that  with 
increasing  number  of  expansion  terms,  the  resulting  loci  of  either  condition 
approach  the  proper  one  more  closely.  The  additional  concentric  hyperbolic 
surface  loci  of  |e  ^  x  E*,|  =  0  exterior  to  the  proper  surface  locus  are  again 
spaced  equidistantly  with  kAR  =  7r/2  in  the  specular  direction,  and  no 
hyperbolic  surfaces  were  found  within  the  proper  one,  though  for  the  larger 
prolate  spheroid  such  surfaces  should  be  expected  but  the  associated  minima  may 
be  too  small  to  be  found  with  the  employed  searching  increment. 

Considering  the  results  exterior  to  the  conical  section  about  specular 
direction,  those  deviations  have  been  expected  by  theory.  It  was  shown  (7644-1-F, 
8579- 1-Q)  that  the  spherical  partial  wave  representation  of  the  near  scattered 
field  was  convergent  down  to  the  minimum  sphere  enclosing  the  loci  of  the  prolate 
spheroid.  Hence  to  determine  the  flank- parts  (equatorial  cross-section)  of 
the  prolate  spheroid,  either  an  alternative  representation  must  be  used,  or 
else  in  the  partial  wave  expansion  the  origin  of  the  coordinate  system  must  be 
translated  from  the  center.  However  in  the  representation  employed,  where  the 
origin  was  at  the  center,  the  results  of  the  condition  f  E  .  -  E  J  1  for  the 

case  ka  =  1,  a/b  =  2  yield  very  satisfactory  results  in  this  region 
(see  Figs.  4-4f  and  4-4g)  if  the  local  polarization  properties  are  such  the  E 
is  parallel  to  the  scattering  surface  in  the  equatorial  plane.  This  behavior,  namely 
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0  =  90° 


FIG.  4-5g:  RESULTING  SURFACE  LOCI  FOR  THE  PROLATE 

SPHEROIDAL  TEST  CASE  ka  =  10,  a/b  =  5,  N  =  16.  (0  =  90°) 

( - ‘{l^il  "  |Eg|}  - {"|EtxE^|2J - {proper  locus 
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12,  16.  0  =  0°,  9  =  135°  . 
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FIG.  4-7d:  PLOT  Of{||E  t|  -  |ES||)  VERSUS  RADIANT  VECTOR  X  =  kR 
FOR  THE  PROLATE  SPHEROIDAL  TEST  CASE  ka  =  10, 
a/b  =  5  FOR  A  VAYING  NUMBER  OF  EXPANSION  TERMS 
n  =  6,  10, 12,  16.  0  =  0°,  9  =  1350  . 
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that  both  boundary  conditions  will  yield  the  best  result  for  these  particular 
polarization  properties,  was  expected  (8579-3-Q)  and  was  found  for  all  employed 
test  cases.  Thus  the  simultaneous  application  of  both  conditions  at  low  frequencies, 
such  that  ka^l  ,  may  be  sufficient  for  a  perfectly  conducting  scatterer  of 
prolate  spheroidal  shape.  The  additional  minima  of  •  E  x  E3^  r  that  lie 

t  * 

within  the  proper  surface  locus  are  of  extremely  small  amplitude  with  respect 
to  that  about  the  proper  surface  locus  and  show  a  sharp  dip,  thus  can  be  determined 
only  if  an  extremely  small  searching  increment  is  used,  and  may  therefore  be 
disregarded. 

Summarizing  the  results  obtained  for  the  prolate  spheroidal  test  cases, 
it  can  be  concluded  as  well  that  the  simultaneous  application  of  the  two  conditions 
at  two  different  operational  frequencies  is  the  proper  method  of  determining 
the  point  which  may  lie  on  the  surface  of  the  scatterer  within  that  conical 
section  which  is  exterior  to  the  minimum  sphere  enclosing  the  two  focal  points. 
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